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Abstract 

We exhibit a static, cylindrically symmetric, exact solution to the Euler-Heisenberg field 
equations (EHFE) and prove that its effective geometry contains (optical) black holes. It 
is conjectured that there are also soliton solutions to the EHFE which contain black hole 
geometries. 
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Chapter 1 

Introduction and overview 



The underlying motivation of the present thesis is the idea of a mathematical connection 
between solitons and black holes. Connections of this sort have been considered elsewhere, 
though not in the same context as the present thesis [1] . We will study a particular system 
of nonlinear PDEs, arising from the Euler-Heisenberg field equations (EHFE), which we 
conjecture has solutions uniting solitonic and black hole-like properties. The EHFE derive 
from the Euler-Heisenberg effective Lagrangian for quantum electrodynamics; the solutions 
that we are interested in are optical black holes. 

1.1 Solitons 

The theory of solitons arises from the study of wave phenomena in nonlinear PDEs. A soliton 
is a solitary traveling wave that maintains its shape through time. Of particular importance 
for us is the nonlinear Schrodinger equation (NSE). In 1-dimension, for a complex wave 
amplitude ip{t,x) with coupling constant 7, the NSE has the canonical form (e.g., Sulem 
and Sulem [2J pp. 5, 20, Drazin and Johnson [3J pp. 34 - 35): 

idtij + dlij + -f\tlj\^ij = 0. (1.1) 

1.2 General relativity 

A spacetime is a 4-dimensional pseudo-Riemannian manifold. The key equation in general 
relativity is the Einstein field equation (EFE), which unfolds to give a system of nonlinear 
PDEs. Solving these PDEs allows one to express the metric coefficients Qf^i, of spacetime in 
terms of the stress-energy tensor T^^. With the cosmological constant A included, the EFE 
reads (Hawking and Ellis [1] p. 74): 

G^u + Ag^u = SnQT^,, (1.2) 

where Q is Newton's gravitational constant and the speed of light is set equal to 1. The 
Einstein tensor G^i, can be expressed in terms of the Ricci tensor R^^, the scalar R = R^, 
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and the metric gf^^: 



R 

Gf^u = R^u '^Qfif (-'-•3) 



Among the known exact solutions to (1.2) are those describing black holes. 

Definition 1. A black hole is a region of spacetime where future- directed outgoing null 
iesics cannot escape. 



Definition 2. A white hole is a region of spacetime where future- directed ingoing null 
geodesies cannot enter. 

Definition 3. The boundary of a black (or white) hole is called an event horizon. 

We use these definitions even outside the context of general relativity. For us, any 
pseudo-Riemannian manifold will be called a spacetime whether it satisfies the EFE or not, 
and one can ask whether or not a given spacetime has black holes. 



1.3 The Euler-Heisenberg field equations 

Quantum electrodynamics can be approximated as an effective field theory governed by the 
Euler-Heisenberg Lagrangian (cf. Euler and Heisenberg [5], Schwinger [6], Novello [7] p. 
292, Boer and van Holten 0): 

L = --F + —(f^ + -gA, (1.4) 
4 90 V 4 y ' ^ ^ 

where a is the fine structure constant, and F and G are the Poincare invariants of the 
electromagnetic field. The speed of light, the reduced Planck constant, the mass of the 
electron, and the permittivity of free space are here set equal to 1. Applying the principle 



of least action to (1.4) leads to the Euler-Heisenberg field equation (EHFE): 

V^F'^^ = — V/, (4^^^^^ + 7GF*'^^) . (1.5) 

Here, F^^ is the electromagnetic field tensor, F*^'^ is its dual, and represents the covari- 
ant derivative with respect to the coordinate x^, using the connection determined by the 



background spacetime metric. Note that (1.5) is, in general, yet another system of nonlinear 
PDEs. 

According to the Euler-Heisenberg effective field theory, the vacuum behaves like a non- 
linear physical medium. Light rays passing through electromagnetic fields are bent as if they 
were passing through water, thus affecting the apparent geometry of objects. This motivates 
the idea that the effective field theory can be interpreted geometrically. Indeed, Novello [7] 
has shown in a seminal work that light rays (small disturbances traveling through the field) 
follow null geodesies with respect to a spacetime metric 'g^y distinct from the background 
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metric ^f^j^. This is called the effective metric. The inverse (or cometric) g'^" can be expressed 
in terms of the stress-energy tensor of the electromagnetic field: 

gi^u ^ Ag^^ + BT"', (1.6) 
where A and B are special functions of the Poincare invariants F and G (see Equation 



(3.57)). This geometrical interpretation of effective field theory demonstrates an analogy 



between nonlinear optics and general relativity, with Equation (1.6) playing the role of the 



Einstein field equation (1.2) [?]. Let us note two subtleties. (1) The effective metric is 
uniquely determined only up to a conformal factor. (2) Since light in the nonlinear vacuum 
experiences birefringence, a given electromagnetic field actually carries two distinct effective 
metrics; one for each polarization state. 

1.4 The idea 



Our main proposal is that the EHFE ( 1.5 ) has soliton solutions with a corresponding effective 



geometry containing a black hole. In this sense, the EHFE would be somewhat in between 



the NSE (1.1) and the EFE (1.2). We have a theorem and a conjecture: 



Theorem 1. There is an exact static solution to the Euler-Heisenberg field equations where 
the effective geometries of each polarization state have black holes. 

Conjecture 1. There is an exact solution to the Euler-Heisenberg field equations which is 
a soliton and whose effective geometries have black holes. 



The theorem is proven in Section |5.3.1[ The soliton of the conjecture is, we believe, an 
imploding solitonic wave. Evidence for this belief is discussed in the section below. The 
purist will note that the hypothetical "imploding soliton" cannot strictly be a soliton as a 
soliton does not change its shape. As the soliton in the conjecture implodes, it will become 
more concentrated and lose its initial shape. However, in cylindrical coordinates {t,r,6,z), 
for a cylindrically symmetric wave approaching the axis r = 0, we conjecture that its radial 
cross section will keep its shape if multiplied by r. 

1.5 Evidence for the conjecture 

Evidence in support of the conjecture comes from nonlinear optics [9]. In particular, Soljacic 
and Segev [lOj have examined the behavior of a beam resulting from the perpendicular 
collision of two plane waves in the Euler-Heisenberg vacuum. Using approximations, they 



determined that the amplitude of the resulting beam satisfies the NSE (1.1). Since an 
imploding wave can be thought of as a limiting case of infinitely many colliding plane 
waves, it seems reasonable to expect that the NSE should be obtained in the case of an 
imploding wave. 

In the work of Brodin et al. [TT], which nicely complements Soljacic and Segev's paper, 
a beam guided between two parallel conducting planes is studied. It was found that the 
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amplitude of this beam satisfies a 2-dimensional cylindrically symmetric NSE. According to 
Brodin et al. [11], for a beam with a certain critical intensity /c, the dispersive and self- 
focusing effects exactly balance and the beam forms an optical soliton of constant width. If 
the intensity I of the beam is less than /c, then the beam width diffracts without bound. If 
/ > Jc, then the beam width collapses to zero in a finite time. These results from nonlinear 
optics show that there is an authentic link between the EHFE and the NSE, which at least 
partly supports Conjecture [1} 



Another piece of evidence for the conjecture comes from the work of Section 5.2 The 



Maxwellian approximation, although it is only a first-order approximation to a solution 
to the EHFE, it gives information on the coordinate velocities of effective geodesies up to 
second-order (see Theorem [s]). When we look at the coordinate velocities of the outgoing 
geodesies to second-order, we find that they are trapped within a certain radius. (There is 
a black hole.) 

Since Conjecture [T] concerns solutions of a nonlinear variational problem, we suspect that 
its proof will use tools from Morse theory (i.e., the calculus of variations in the large). 



1.6 Organization 

This thesis is organized as follows. 

Chapter|2]is a pedestrian introduction to the required mathematical physics. Chapter|3]is 
a self-contained review of Novello's theory of effective geometry. In Chapter |4] (which can be 
omitted on a first reading), we study the effective geometry of plane waves, and calculate the 
index of refraction through a plane wave confirming earlier approximations done by others 
using different methods. In Chapter [5} we use well-known solutions from Maxwell's theory 
(for imploding cylindrically symmetric waves) and investigate the corresponding effective 
geometries which resemble black holes. At the end of Chapter |5j we prove Theorem [T] by 
explicitly finding an exact solution to the EHFE with the required properties. Although this 
exact solution is static, it shares some similarities with an ingoing cylindrical wave solution 
because its Poynting vector points radially inward. 
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Chapter 2 
Preliminaries 



The present chapter is meant to be a self-contained pedestrian introduction to the relevant 
mathematical physics. 



2.1 Nonlinearity of the vacuum 

According to quantum electrodynamics, photons can scatter off of each other. This photon- 
photon scattering effect, also known as the nonlinearity of the vacuum, was calculated by 
Euler and Heisenberg [S] in the mid-1930s, but it is so subtle that no currently available 
experiment is yet sensitive enough to measure it. 

Photon-photon scattering arises from processes that involve virtual electron-positron 



pairs (see Figure 2.1). In the classical limit, the effect of these virtual particles on real 
photons can be approximated by introducing nonlinear terms to the Maxwellian Lagrangian. 
In this so-called effective field theory, the photons do not necessarily follow null geodesies 
in the background metric. Instead, they follow null geodesies with respect to a so-called 
effective metric, as will be explained in Chapter |3| 



2.2 The background spacetime 

In a curved background, according to Drummond and Hathrell [12], the physical Lagrangian 
for the effective field theory acquires a nontrivial dependence on the spacetime curvature. 
In the present thesis, we restrict ourselves to a Minkowskian background, so these curvature 
coupling effects can be ignored. 

The background metric tensor is denoted by g^^. Its inverse, the so-called cometric, 
is denoted g'^'^ and one has that g^^gxi, = 5^^, where 5^^ is the Kronecker delta (and the 

Einstein summation convention is followed as usual). The present work uses the H 

signature convention. 

Recall that a given metric determines a unique torsion-free connection V by requiring 
that the covariant derivative of the metric be zero (e.g. Hawking and Ellis p. 40, or 
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Figure 2.1: Feynman diagram of photon-photon scattering. 



Spivak [13J pp. 236 - 237): 



^\9ixu = 0. 



(2.1) 



The connection coefficients (Christoffel symbols) tlius determined are given by: 



where the operator denotes partial differentiation with respect to the x'^ coordinate. I.e., 
■= d/dx'^. 

2.3 Electromagnetic fields 

The electromagnetic field F^^, is a closed 2-form, and locally there exists an A-field such 
that: 



(2.2) 



F 



d^A, - d,A^ (= V^A, - V, 




(2.3) 
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Note the Bianchi identity: 

= dxF^, + 9,Fa^ + d^F^x 



(2.4) 



The antisymmetrization of an indexed quantity is indicated by placing square brackets 



around indices. For example, given the tensor T^^ 



we can write: 



(2.5) 



where the sum is taken over all permutations a of the indices ai, a2, ■ ■ ■ , a^. The value of 
sgn((T) is +1 if cr is an even permutation of the sequence a\a2 ■ ■ ■ O-n and is —1 if a is an odd 
permutation. 



Using antisymmetrization, Equation (2.4) can be written as: 

= 9[aF^.] (= V[aF^.]). 
The Levi-Civita tensor e^^^y is defined such that: 



(2.6) 



(2.7) 



where g denotes the determinant of the metric g^y. Note that the value of Sa^^y is 
unless the indices a,/3,/i, are all distinct. Furthermore, note that Ea^y^y is +\/|^ if the 
sequence a/S/j^iy is an even permutation of the sequence 0123 and is — a/M if oif^Hi^ is an odd 
permutation. 

The Levi-Civita tensor allows us to express the Hodge dual of F^y by writing: 



{21 



Unless otherwise specified (e.g., in Section 3.3), indices are always raised or lowered with 
respect to the background metric. So, e.g., F^^ = g'^^g^^ Fas- 



Note that the Bianchi identity \2A\ can be expressed in terms of the dual tensor by 

(2.9) 



writing (cf. Landau and Lifshitz [H] p. 67): 



= d^.F*^" 



V^F*^'') 



2.4 Effective Lagrangians 

The physical behavior of the electromagnetic field is governed by a Lagrangian L which is a 
scalar function of the field A^, its covariant derivatives Vfj,Ay, and the background metric. 
The field equations are obtained from the principle of least action: 



dL _ dL 

'diV^Ay) ~ My- 



V^ ^,^ ^ , = — . (2.10) 
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In the case of Maxwell's theory, the Lagrangian in the absence of charges and currents 
(a field in the vacuum) can be written out as: 



[V^A, - V.A^) {V^Ap - VpA^). (2.11) 



If we define F := F^^F'^", then Equation (2.11) simplifies to L = —F/A. Using (2.10), one 



recovers the familiar form of Maxwell's equation for the vacuum: 

Vf.F^" = 0. (2.12) 

The quantity F introduced here is a scalar invariant of the electromagnetic field tensor. 
In fact, there are only two algebraically independent scalar invariants for F^^. These are 
represented by the so-called Poincare invariants F := F^i.F'^'^ and G := F'^'^F*^ (Landau 
and Lifshitz [II] p. 64). 

We define the class of L{F,G) -theories as electromagnetic theories in which the La- 
grangian L can be expressed in terms of the Poincare scalars F and G alone, i.e., L = 
L{F, G). It is assumed that the partial derivatives of L = L{F, G), with respect to F and G, 
exist at least up to second-order, and that they are continuous. We will use the notations 
Lp := dL/dF, Lg := dL/dG, LpF ■= d^L/dF^, Lgg ■= d^L/dG^, Lpc ■= d^L/{dGdF), 
etc. We observe that L{F, G)-theories are guaranteed to be Lorentz invariant since both 
F and G are Lorentz invariant quantities. A particularly important L{F, G)-theory is the 



Euler-Heisenberg theory (to be introduced in Section 2.5). 

Chapter [3] is primarily concerned not with a particular theory but with general L(F, G)- 
theories. However, in subsequent chapters, attention is restricted to the Euler-Heisenberg 
theory. 

Although effective field theories more elaborate than the L{F, G)-type can be constructed 
by writing Lagrangians that include terms involving the covariant derivatives of the field 
(e.g., terms like V\F''^^V^F^'^) we will not deal with such things in the present work. 



2.5 Euler-Heisenberg theory 

Euler-Heisenberg theory is an effective field theory which approximates the physical theory 
of quantum electrodynamics in Minkowski spacetime. 

Up to second order in the fine-structure constant a, the Euler-Heisenberg Lagrangian 
is given by (cf. Euler and Heisenberg [5J, Schwinger [6J, Novello (7j p. 292, Boer and van 
Holten [H]): 

L = -^F + —(f^ + -gA. (2.13) 

4 90 V 4 y ^ ^ 

Since we are presently working in natural units, the speed of light c, the reduced Planck 
constant h, the mass of the electron rrie, and the permittivity of free space eo, are here set 
equal to 1. Conventional Lorentz-Heaviside units can be restored by replacing with 
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a'^h^'m~'^c~^. (Sometimes we use the letter a as a tensor or pseudotensor index, but no 
confusion between a as the fine-structure constant and a as an index should arise because 
the context will make the meaning of a clear.) 



Equation (2.13) applies to fields having strength A and frequency oo such that [S]: 



A < 

v47ra 

w < 1. 



(2.14) 
(2.15) 



In other words, the field strength should be much weaker than the critical field 1/ a/ Aira and 
it should be approximately constant on scales much less than the Compton wavelength of 
the electron (which is unity, in our units). In the present work however, we will not worry 



about these physical restrictions (2.14) and (2.15). 



2.6 Field equations 



Theorem 2. For a Lagrangian of type L = L{F,G), the principle of least action (2.10) 
yields the field equations: 

{LfF^"" + LgF*^") = 0. (2.16) 



Proof. For a Lagrangian of type L = L{F,G), we get that dL/dAy = 0. So Equation (2.10) 
reduces to 

dL 



We get that: 



because: 



and: 



dL ^ OF ^ dG , , 

Lf ^,^ , , (2.18) 



d{W^A,) 'd{V,A,) "5(V^A,)' 



dF 

AFf"", (2.19) 



d{V^A, 



AF*^\ (2.20) 



d(y^A,) 

When an index is free and not to be summed over by the Einstein convention, we draw a 
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bar over it: 



dF 



d 



d 
d 
d 



[g'-^g^'F^pF^,) 

{g"'g^' {V^Ap - VpA^) {VxAp - VpA^)) 

{g'^'g'' (V^Ap - v,Ap) {VxA, - VpA^) 
+ g'^g'^' ( VpA^ - VpA,) {VxAp - v,A^) 



-g'^g'^" (VaA,-VpA) 

jpp.P _ pup. _|_ pfiU _ pup. 

iF'^''. 



VpA^) {VpAp - V^Ap) 



Equation (2.20) can be established similarly. 



The result (2.16) follows from Equations (2.18) - (2.20). 



□ 



Specializing (2.16 ) to the Euler-Heisenberg Lagrangian (2.13 ), we get the Euler-Heisenberg 
field equations: 



a 



V. (AFF^" + 7GF 
45 



*fiu\ 



(2.21) 



For an arbitrary G')-theory, Equation (2.16) implies: 



= V^iLpFn + '^ALaF^n 

= {LFFV,F + LFGy,G)F^- + LFV^F^^^ 



(2.22) 



The last term is zero by the Bianchi identity (2.9). By computation, one notes that V^F 
2F°^V/,F„;3 and V^G = 2F*"^V^F„/3. Hence, if we define the tensor [15J: 



Q 



afifiu 



LffF^^F^"" + Lfg (f^'^F*"'' + + LggF*''^ F*^"" , 



(2.23) 



then we can rewrite the field equations (2.16) as: 



0. 



(2.24) 



Assuming Lp ^ 0, Equation (2.24) can be rearranged: 

2 



(2.25) 



(The case Lp = Q is discarded since it is not physically interesting.) 

2.7 Stress-energy tensor 

Given a Lagrangian L, one can define a stress-energy tensor T^y through the equation: 

dL 



T ■= 2 



(2.26) 



This expression for the stress-energy tensor is implicit in e.g. Novello [7J pp. 271, 275, 
Landau and Lifshitz p. 77, Hawking and Ellis jl] p. 66, and Poisson [16] p. 125. 
Observe that different authors disagree on the overall sign on T^^ due to differing signature 
conventions for the metric. 



For a Lagrangian of the form L = L{F, G), Equation (2.26) gives: 

dF 



T 



2 L 



7 



+ L 



dG 



G 



dg 



-Lg 



—ALt?F'^F —ALr^F'^F* 



Lg, 



(2.27) 



Using the well-known identity 4F^"F*^ 
agreement with Novello [7j p. 275): 



-Gg^^ (cf. Novello [7] p. 272), we get (in 
T^, = -4L^F/F„, - (L - GLG)g,u. (2.28) 
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Chapter 3 

Effective geometries 



The purpose of this chapter is to give a quick self-contained review of Novello's theory 
of effective geometries in nonhnear electrodynamics [3 [151 113 HI [IS]- Our exposition is 
informed by the existing literature, most notably the work of Novello [7]. However, we do 
not follow any specific work too closely. 

3.1 Electromagnetic shock waves 

The wave front of an electromagnetic shock wave is defined by a hypersurface S across 
which the field derivatives are discontinuous. Given a set of local coordinates x'^ for the 
(background) spacetime manifold, this hypersurface S can be described as the set of solutions 
to the equation: 



We will need to assume that the first-order partial derivatives of z{x^) exist and are con- 
tinuous on S, and that the gradient /c^ := (9^2 does not vanish on S. The hypersurface 
S, at least locally, splits the manifold into two regions := {x^ : > 0} and 

dJl- := {x^ : ;z(x^) < 0}. 

The jump of an arbitrary function J through S is denoted by the Hadamard bracket 
[J]e. For each point p of S, we define: 



where the points p"*" and which tend towards p, belong to the regions OJt"'" and 
respectively (Papapetrou [20] p. 170). Note that if J is continuous across S, then [J]s = 0. 
The converse is not strictly true, since it is possible to have a function with a so-called 
simple discontinuity whereby \mip+^p J{j)^) = limp-_^p J(p^) 7^ J{p)- On the other hand, 
the derivative of a function can be discontinuous but the discontinuity is never of the simple 
type (Rudin [21] p. 109). In light of this, a partial derivative d^J is discontinuous across E 
if and only if [<9^J]s 7^ 0. 



zix^") = 0. 



(3.1) 



[J]s(p) := lim J{p'^) — lim J{p ), 



(3.2) 
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Since E is the front of an electromagnetic shock wave, the electromagnetic field is con- 
tinuous across S but some of its derivatives are discontinuous across S. We express this by 
writing: 



lj.u\T, 



0, 



and: 



\^xF^y\j] 7^ for some A, /i, v. 



(3.3) 



(3.4) 



Similar conditions hold for the dual tensor F, 



Note that since the field and the Christoffel symbols T^^^ are continuous, we have 



Now consider a second coordinate system {x^} such that 



[20] p. 171). Then S can be described by the equation x° 
continuous first-order partial derivatives, we get that: 



mHs- (3.5) 

z{x^) (cf. Papapetrou 
0. Moreover, since zix'^) has 



dx^ 

[{dQF^^)dxz]^ 
[<9oi^/.i.]s ■ dxz 



(3.6) 



where /^j^ := [^gF^i^jE is the so-called discontinuity or disturbance in the field, and the 
1-form kx '■= dxz is called the propagation vector. It is required that kx be nonzero. 

Theorem 3. The quantity f^y is a tensor. Moreover, it is a 2-form. 

Proof. We verify that if we go to another coordinate system {x^'}, then the quantity f^i, 
transforms tensor should. 



ffi'. 



do 



fit/ 



dx^ dx^ 
dx^^' dx^' 



dx^ dx" ( dx^ dx^ 

7 + -r/x!/0'o 



^x^^' dx'- 

dx^ dx" 
dx^^' dx"' 

/ dx^ dx" 
^ \^x^'' dx"' 
dx^ dx" 



dxi^' dx"' 



J s 



+ 



J s 



. dx^ dx" 
^ dxf^' dx"' 



continuous 



dx"' 



(3.7) 



21 



Moreover, f^i, is a 2-form since f^i, = —fu^i- (Note that the above underhned "continuous" 
term is continuous since the background spacetime, which is Minkowskian, is C^.) □ 



Note that for the dual F*^^ we write, in analogy with Equation (|3.6|): 



[dxF;,h = f;ukx, (3.8) 

where f*^ := [(9o'-F*j,]s is the discontinuity of the dual field. Analogously to the relation 

Kfi = I^ais/.^F'"', one has: 

= \eo.p,ur. (3.9) 

Moreover: 

[dxFn^ = rh, (3.10) 

and: 

[SaF^He = r^^kx. (3.11) 



3.2 Dispersion laws and polarization 

In nonlinear field theory, field discontinuities (or photons, in a classical corpuscular sense) 
can exhibit birefringent behavior [7l[T7l[T5]. This means that the way a photon propagates 
through the field depends on its polarization state. Whether a theory predicts birefringence 
or not depends on the Lagrangrian used. For example, in Born-Infeld electrodynamics, there 
is no birefringence (see e.g. Novello [7] p. 276). In the Euler-Heisenberg theory however, 
there is. 

The goal of the present section is to derive the dispersion laws for L(F, G')-theories. We 
begin with the following observation: 

Theorem 4. Locally, there exists a l-form such that /^j^ = Pf^k^^ — p^k^. 



Proof. Applying the Hadamard bracket to both sides of Equation (2.9) gives: 



r^k, = 0, (3.12) 

which implies that det (/*'''') = -\g\-^/^ dei{f^^) = 0. Since we are in four dimensions, it 
follows that f^u is simple (e.g., Penrose and Rindler [22] p. 166). That is, locally there exist 
1-forms and such that: 

f^^u = U[^v^] = ]r{u^v^ - u^v^). (3.13) 



Taking the Hadamard bracket of the Bianchi identity (2.4) gives: 



fl^ukx] = 0. (3.14) 
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Equations (3.13) and (3.14) imply that the triple wedge product of u^, Vy and k\ vanishes. 



Hence m^, Vy and kx must be linearly dependent (e.g., Madsen and Tornehave [23] pp. 11 
12) and so locally there exists a 1-form for which we have the decomposition: 



(3.15) 

□ 



Theorem |4|/Equation (3.15) says that the field discontinuity f^i, (or photon, as we are 
apt to call it) is the wedge product of the propagation vector k^ together with p^. Without 
loss of generality we can assume that is orthogonal to k^ and thereby interpret as 
being the (non-normalized) polarization vector (actually a 1-form) for the photon. 



Let us take the Hadamard bracket of both sides of Equation (2.24). After a bit of 
rearranging, one gets that: 



(3.16) 



Substituting (3.15) into (3.16), and using the assumption that p^ is orthogonal to k^, it 
follows that: 



Lpk^P^ = -^QT^k^kpp,, 
where /c^ := g'^'^k^ku- Assuming Lp 7^ 0, we can write: 



(3.17) 



'Qri ^kfykpPp 



For convenience, define the tensor (cf. De Lorenci et al. [T5] ' 



(3.18) 



Then Equation (3.18) can be expressed as: 

S^P, = 0. 



(3.19) 



(3.20) 



If k"^ 7^ 0, Equation (3.18) implies that p^ can be expressed as a linear combination of 
:= F^^k\ and h*^ := F* kx. (Note that both /i^ and h*^ are orthogonal to k^: since F^'^ 
and F*'^" are skew-symmetric we get that h^^k^ = F^^kxk^j. = and h*^k^ = F*'^''kxk^ = 0.) 



Writing: 



We get that: 



(3.21) 



4 



Lpk'^ 



+ LFpe + Lpch^'h*^ K + {LcGh'^K + Lpch'') h*^ , (3.22) 
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and 



Lpk^ 



+ LGGh*''hl + LpGh'-hl] hi) .{3.23) 



the well-known identities (cf. Novello [7] p. 272): 



Equations (|3.22|) and (|3.23|) can be recast into a somewhat more useful form by exploiting 

(3.24) 



and: 



For (3.24) and (3.25) respectively, contracting both sides with k^ky gives: 

1 



and: 



h^hl = -Gk\ 



h*"hl + h'^ = -FP. 



Using (3.26) and (3.27), Equations (3.22) and (3.23) become: 



4 



+ -GLfg I A;' + Lpph'' \K+[ -GLGck^ + Lpch^ I h 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



and: 



4 

+ 



^GLpp - ^FLfg ] + LFGh"" ] K 



e + LGGh' 1 h: 



(3.29) 



Equations (3.20), (3.21), (3.28), and (3.29) give: 



L, 



4 4 



+ -GLfg k"^ + Lpph'' 



^GLff - ^FLfg ] k' + LFGh' ] ] K 



+ ( a ( -GLGGk' + LFGh' 1 + h 



^ - ]^FLgg + \gLfg ) k^ + Lgg^' 1 1 K. 



(3.30) 



24 



First we consider the case where hi, and are hnearly independent. In this case, we 
have the following linear system in the variables a and b: 

a{{^ + IGLfg) + LpFh') + b {{iGLpp - \FLfg) + LpGh') = 

(3.31) 

a [iGLcGk^ + Lpch'') +h{{^- \FLgg + \GLfg) k^ + Lech'') = 0. 

The determinant of this system has to be zero (there is a nontrivial solution for a and b 
because the polarization vector = ah^ + bh*^ is nonzero). Thus: 

((^ + \GL,a^ e + L^,h'^ ((^ - \FLac + \GL,a^ + Lach'^ 

-GLff - ]^FLf(^ e + Lpch^^ {^^GLcGk^ + Lpch''^ . (3.32) 

In the case where /i^ and h*^ are linearly dependent, it follows that S^^^h^ = S^^^h* ^ = 
since = ah^ + bh* ^ and S^^^p^ = 0. We thereby obtain the system: 




((^ + IGLfg) k^ + Lpph^) K + {\GLGGk^ + LpGh'') K = 

{{\GLff - \FLpg) e + LpGh^) K + ((^ - \FLgg + \GLfg) k^ + LcGh'') K = 0- 



(3.33) 



Since we require at least one component of or h* to be nonzero (p^ = ah^ + bh* is 



nonzero). Equation (3.32) holds even if /i^ and /i* are linearly dependent. 

Expanding the products and combining like terms. Equation (3.32) can be put in the 
form: 



where we define: 



Ai 

A2 

As 



{Lp + GLpcY — LggC^FLp + G'^Lpp), 
4 [Lp{Lpp + Lgg) + '^F{LjpQ — LppLgg)) 
16{LppLgg — L 



FG) 



(3.34) 

(3.35) 
(3.36) 
(3.37) 



We now have the following result: 
Theorem 5. Assuming that Li? 7^ 0, A;^ 7^ 0, Ai 7^ and Ag — 4A1A3 > (^m order to 



ensure that Equation (3.34) gives real solutions for k"^), we have the dispersion law(s) l2^ : 



where: 



A±: 



e = k±h\ 



-A2 ± ^{A2? - 4A1A3 
2Ai 



(3.38) 



(3.39) 
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Next, we will show that Equation (3.38 



) continues to hold even if k'^ = 0. More precisely: 



Theorem 6. If k'^ = and LffLgg — Lpc 0? then h"^ = 0. 

(Note that with the Euler-Heisenberg Lagrangian (2.13), we have LffLgg 



^FG 7^ 



Proof. For an indirect proof, suppose that k"^ = and h"^ ^ 0. Note that when k"^ = 0, 



Equations (3.26) and (3.27) give h'^h*^ = and = 



Note that /i^ cannot be timelike 



because otherwise /i* would also be timelike and one cannot have two orthogonal timelike 
vectors in Minkowski spacetime. The only remaining possibility is that /i^ (and consequently 
h* ) is spacelike. We show that this implies LffLgg ~ L%q = 0. 



To this end, note that with k^ = 0, Equation (3.17) becomes: 







Q^'^'^k^kfsp, 

LFph^ (h^Pu) + Lfg (h^ (/i* V) + h*^ {Wp,)) + LGGh*^ ih*''p. 



(3.40) 



Contracting (3.40) with and /i*^ respectively, and using the relations h*"^h*^ = h 



and h'^h*^ 



0, it follows that: 



Lff (h^p^ 
Lfg (h^P. 



+ Lfg (/i* V) = 
+ Lgg {h*'Vu) = 0. 



(3.41) 



We claim that h'^py and h^^p^ cannot simultaneously be zero. To establish this claim, note 
that since k^ is null, and since p^ is orthogonal (and not parallel) to /c^, it must be that p^ is 
spacelike. So if h^Pu and h*^py were both simultaneously zero, we would have an orthogonal 



basis consisting of three spacelike 1-forms h 
possible. 



h* 



and a null 1-form /c^, which is not 



Consequently, the system (3.41) has a nontrivial solution for h^pu and h*'^p,^, and so it 
must have a vanishing determinant: LffLgg ~ L'j^Q = 0. □ 

Now, according to Theorem [s} there can be two possible values of /c^. This has to do 
with the fact that a given photon (field disturbance) is in one of two polarization states. To 
understand why this is true, note that (as explained in e.g., De Lorenci et al. |1^.) for each 
possible value of fc^ there corresponds a certain solution space for the unknowns a and h in 



the system (3.31), and p^ = ah^ + hh* . (Here we identify the solution space of a and h with 



the "polarization state." 
By defining: 



± 



ALff + {Lf + GLfg)^± 
4:Lfg + GLgg^± 



(3.42) 



Equation (3.38) becomes 



jFF 



+ Lfg^± 



Lf + {Lfg + Lgg^±)G 



(3.43) 
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which matches Equation (24) in De Lorenci et al. |15j . 



We note that with the Euler-Heisenberg Lagrangian (2.13), Equation (3.39) reads: 

224«2 



A 



± 



495 + 12Fa2 ^ ^18225 - 18360Fa2 + 4624^2^4 + siseG^a^' 



(3.44) 



3.3 Effective null geodesies 



The dispersion laws in nonhnear electrodynamics have an appealing geometric interpreta- 
tion, where they are thought of as light cone conditions in a so-called effective geometry. 



Using the fact that h = —F^^F°'^kfj_ky^ we can write Equation (3.38) in the form 



Provided that the symmetric tensor 'g^'^ defined by: 



(3.45) 



(3.46) 



is nonsingular, an effective metric g^i, can be defined such that 'g^^'gxv = 5^,^. We get the 
effective geometry by treating 'g^y as if it were the metric for spacetime. 

We can think of /c^ as being null with respect to the effective metric since: 



r^KK = 0. 



(3.47) 



Moreover, the integral curves of k^ (i.e., photon worldlines) turn out to be geodesies 
with respect to the effective metric: 

Theorem 7. The integral curves of k^ are null geodesies with respect to the effective metric. 
Proof. The proof is given in Novello ^ pp. 273 - 274. We repeat it in order to be self- 



contained. The first step is to take the partial derivative of Equation (3.47) to get: 



2{dy,k^)k,g^'' + k^,k,d,:g 



'/ii/ 



0. 



(3.48) 



Next, exploit the fact that the effective metric g'^", provided it is nonsingular, determines a 
set of torsion-free connection coefficients (through the usual Christoffel formulae, i.e.. 



Equation (2.2)) and thereby determines a covariant differential operator Va such that: 



Contracting Equation (3.49) with kf^k,^ one gets: 

kf.kydxg^-' = -2k^kyf\,r''- 



(3.49) 



(3.50) 



By substituting (|3.50|) into Equation (3.48), it follows that: 

"^^^ ^Vxk,)k, := r''(dxk,-f\^k^)k. 



9' 



0. 



(3.51) 
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Since := dfj_z, and since partial derivatives commute, one gets that: 



Defining k^^ := g^'^k^, and using (3.52), Equation (3.51) can be rewritten as: 



which imphes that we have a geodesic (e.g., Poisson fTU] p. 61) 



(3.52) 

(3.53) 

□ 



Note that a given effective metric is only defined up to conformal equivalence because 
all that we have are null geodesies. More precisely, an effective geometry is an equivalence 
class of conformally equivalent Lorentzian metrics. However, we will only work with one 
representative at a time, choosing whichever conformal factor suits our fancy. 



In the case of Euler-Heisenberg theory, the effective cometric (|3.46|) becomes: 

9'^ + 



9' 



495 + 12Fa2 ^ ^18225 - 18360Fa2 + 4624^2^4 + 3136G^a^' 



(3.54) 



3.4 The relationship between the effective metric and 
the stress-energy tensor 

For an electromagnetic field governed by a Lagrangian of the type L = L{F,G), we found 



that the stress-energy tensor is (2.28): 

T„„ = —ALpF^'^Fau — {L — GLg)9^u- 



(3.55) 



Raising the indices of Equation (3.55), and rearranging, we get that the substress tensor 
pii^pau -g (^assuming Lp ^Q): 



a 



4L? 



(T/^- + (L - GLg)9 



(3.56) 



Using Equation (3.56), we can rewrite Equation (3.46) as 



fit/ ^i: rj-i/iu 



(3.57) 



Equation (3.57) is analogous to Einstein's field equation from general relativity in that it 



relates the effective geometry with the stress-energy tensor. 

We note that since the effective metric has this direct dependence on the stress-energy 
tensor, de Oliveira and Perez Bergliaffa [25] have suggested that the Segre classification 
of the stress-energy tensor yields a simple classification scheme for effective geometries in 
nonlinear electrodynamics. 
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Chapter 4 
Plane waves 



The main purpose of this chapter is to investigate the effective geometry of a circularly 
polarized plane wave. Our findings independently confirm those of a 2002 paper by Denisov 
and Denisova [26]. The effective geometry of a plane wave is conformally equivalent to 
Minkowski spacetime but it is distinguishable from the Minkowski background because the 
effective null geodesies are not necessarily null in the background. 

Moreover, we show that, as viewed from the background coordinates, shock disturbances 
in a circularly polarized monochromatic plane wave field propagate with a directionally 
dependent index of refraction which can be easily computed. Our result for the index of 
refraction confirms, to lowest order, the calculations performed by Affleck [27] in 1988. 
Lorentz invariance is manifestly preserved throughout in our approach. This is in contrast 
to Affleck's non-invariant approximation. 

We close the chapter with a discussion anticipating the possibility of optical black holes 
in vacuum. 



4.1 Effective geometry of null fields 

Recall that a null field is one such that + = 0. In the case of a null field, the 



nonlinear field equations (2.16) reduce exactly to Maxwell's equations in the absence of 
charges and currents, provided that Lq = when + = 0. Note that the Euler- 
Heisenberg Lagrangian (2.13) indeed satisfies this latter condition. Consequently, a null 



field is an exact solution to the Euler-Heisenberg equations if and only if it is an exact 
solution to Maxwell's equations. 

Let us consider the effective geometry corresponding to a general null field. 



In Section 3.2, we found that field disturbances (shock waves) disperse according to: 



1 + 



ALp 



9' 



A- 



0, 



(4.1) 



where K± is given by Equation (3.39). The choice of ± depends on the polarization of the 
disturbance. Accordingly, we will refer to the polarization modes as being either + or — 
modes. 
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Specializing to null fields, Equation (4.1) leads to an effective cometric given by: 



9"^ 



where: 



R 



A± 



± 



F2+G2=0 



(4.2) 



(4.3) 



ALF + A±iGLG-L) 

For the Euler-Heisenberg Lagrangian, we get that P± = (22 ± 6)a;^/45 (cf. De Lorenci et 

ai. [ini). 

Locally, the stress-energy tensor for a null electromagnetic field satisfying the dominant 
energy condition (see e.g.. Hawking and Ellis [1] p. 91) can be expressed as |28j: 

where is a null vector. Thus effective cometrics corresponding to null fields can be 
expressed by equations of the form: 

= gt^^ + Pjn\ (4.5) 



4.2 Effective geometry of plane waves 

Taking the usual t, x, y, z coordinates for Minkowski spacetime. The field tensor: 
/ F'tt Ftx Fty Ftz \ 



HI/ 



F , F F F 

^ xt ^ XX ^ xy ^ xz 

F t F F F 

^ yt ^ yx ^ yy ^ yz 



zt zx 



zy 



Fzz J 



( Aco^{uj{t - z)) 

-Acos(a;(t - z)) 
-5sin(w(t - z)) 



5sin(a;(t - z)) \ 

A cos [bjit — z)) 

5sin(a;(t - z)) 



v 







-A cos (a;(t — z)) — i? sin (a;(t — z)) 







/ 

(4.6) 

describes a monochromatic plane wave propagating along the +2; direction, having frequency 
w, and elliptical polarization with fixed amplitudes A and B. The stress-energy tensor 
corresponding to this field is: 

/ rpti rptx rpty rptz \ 
rpxt r-pxx 'pxy rpxz 



rpflU 



ipyt 'pyx 'pyy 'pyz 
^ rpzt rpzx 'pzy rpzz J 

f A^ cos\oo{t - z)) + sm\ooit - z)) A^ cos\uj{t - z)) + B^ sm\uj{t - z)) \ 





\A^cos^{Loit-z)) + B^sm\Loit-z)) A^ cos^{ujit - z)) + B^ sm\uj{t - z)) J 

l^^l", (4.7) 
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where: 



/ \ 

ly 

\p J 

( ^A^ cos2(a;(t - z)) + sin^{uj{t - z)) \ 




V ^A^ cos2(u;(t - z)) + 52 sm2(w(t - z)) / 



(4i 



Note that in the case of circular polarization {A = B), the stress-energy tensor takes on 
a particularly simple form; it becomes constant: 



rpfiu 



( A^ 








A^ \ 


























A^ 








A^ y 



(4.9) 



Since the stress-energy tensor (4.9 ) is constant, the corresponding effective geometry (4.2 ) 
must be flat. In the case of an arbitrary elliptically polarized plane wave, the stress-energy 
tensor is no longer constant but one can nevertheless calculate that the Riemann curvature 
of the effective geometry still vanishes identically. This confirms results first published 
by Denisov and Denisova in 2002, who found that the effective geometries corresponding 



to monochromatic plane waves are flat using the Euler-Heisenberg Lagrangian (2.13) 
More generally, we note that if we have a Lagrangian of the form L = L{F, G) and if Lq = 
when + = 0, then the resulting field theory will have plane waves as exact solutions 
and these will yield fiat effective geometries. 

Although the effective geometry of the plane wave is fiat, the effective null geodesies 
are not necessarily null with respect to the fiat background metric (the precise manner in 
which the effective light cones embed in the background is studied in Section 4.4). Thus 



we see that effective geometries which are fiat can nevertheless be distinguishable from the 
flat background spacetime. This phenomenon is not limited to plane waves or even to null 
fields. Just to give a concrete example, a constant uniform electric field (e.g., in a rest frame, 
[At, Arc, Ay, A z) = {0,0,0, Et) with E constant), satisfies the Euler-Heisenberg equations 
(2.21) and its effective geometry is a copy of Minkowksi spacetime. Indeed if the stress 



energy tensor of a field is constant then its effective geometry must be flat. Considering the 
symmetry of such a situation, this should be expected. 



4.3 Refraction in a circularly polarized plane wave 

The purpose of this section is to calculate the index of refraction for field disturbances 
(low- intensity external "photons") propagating in a circularly polarized plane wave field. 
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Note that the stress-energy tensor of a null field given by two constant perpendicularly 
crossed electric and magnetic fields is the same as that of the circularly polarized plane wave 



(4.9). Thus, the effective geometry of crossed null fields is the same as that of the circularly 



polarized plane wave. 



Using the stress-energy tensor (4.9), together with Equation (4.2), we calculate that the 



effective metrics corresponding to a circularly polarized plane wave (propagating in the +z 
direction, with amplitude A) are: 



;i - P±A^)dt^ + 2P±A^dtdz - dx^ - dy'^ - (1 + P±A^)dz^ 



(4.10) 



Since the Christoffel symbols for the effective metric (4.10) vanish identically, the null 



geodesies in the effective geometry are simply rectilinear curves in the coordinates t, x, z. 

Consider an effective null geodesic that passes through the origin of the coordinates 
t,x,y,z. The corresponding projection (i.e., light ray) for this geodesic in the three- 
dimensional X, y, z space issues from the origin and intersects the unit sphere at spherical 



coordinates 6 by (p (see Figure 4.1). The angle (p measures the angle that the ray makes 
with respect to the +z measured in the x, y, z system. 



Z 




Figure 4.1: A light ray issuing from the origin intersects the unit sphere at 6 by ip. 
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Using the standard conversion formulae between rectilinear and spherical coordinates 
{x = rsinipcosO, y = rsinyjsin^, z = rcosip), Equation (4.10) implies that, along an 
effective null geodesic through the origin: 



(l + P±A'^ cos^ ip) 



dz\' 



dz\ 
dt) 



Thus: 



dz 
'dt 



2P±A^ cos^ — - (1 - P±A^) cos^ = 0. 



(4.11) 



P±A^ cos^ LP + cosv^Vl - P±A'^ sin^ 



1 + P±A^ cos'^p 



(4.12) 



Consequently, as measured in the t, x, y, z coordinates with respect to the background metric, 
discontinuities in the plane wave field propagate with a y9-dependent velocity: 



v{p>) 



dx'^ + dy"^ + dz'^ 



df^ 



P±A^ cos p> + y/l - P^A^ sin^ p> 



l + P±A^ cos2 p> 
So the plane wave has an index of refraction, n{ip) = l/v{p>): 

1 + P^A'^ cos^ ip 



(4.13) 



n{(p) 



P±A^ cos <p + Vl - P±A2 



sin^ {p 



(4.14) 



Expressing Equation (4.14) as a power series in A, one gets: 

n(<^) = l + 2P±A2sin^(|) +P20(A^ 



(4.15) 



Affleck [27] approximated a formula for n{p>) using methods different from ours. The formula 
which he obtained (correcting for typos) is nothing but the first two nonzero terms in the 
expansion (4.15). (Note that Affleck's formula for n{p>) was apparently published with a 
small typo: in his Equation (14), the factor [eEQ/rn?) should be {eEQ/im?Y .) 



4.4 Visualizing effective light cones 



The purpose of this section is to describe how the light cone structure of the effective 



geometry given by Equation (4.10) embeds in the background spacetime. 



From Equation (4.13) we get that, as seen in the background, the plane wave induces 
a drag effect for field disturbances. Low-intensity photons that probe the field along the 
direction of the plane wave (the direction given by the so-called Poynting vector) will con- 
tinue to travel at the usual speed of light: f (0) = 1. Along any other direction, the field 
disturbances are made to travel at less than the speed of light. This drag effect is most 
pronounced for = vr, the direction exactly opposite to the Poynting vector. 
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According to Equation (4.13), when P±A'^ > 1, field disturbances are confined to prop- 



agate only in directions such that csc"^ ip < P±A'^. 

If a shock wave issues from the origin, then we can calculate the location of the wavefront 
in x,y,z space after a unit t-time by plotting Equation (4.13) in the a:2;-plane using {v,ip)- 
polar coordinates (i.e., z = vcosip and x = vsimp). One can then rotate this graph about 
the z-axis {ip = 0) in order to visualize the wavefront as a surface of revolution (in fact, 
the surfaces in this case are ellipsoids). Plots of (4.13), representing the range of qualitative 



behaviors, are given in Figure [4^ Note that disturbances propagating against the direction 
of the field experience a kind of drag effect. The case A = (solid red in the figure) is a 
limiting case in which the plane wave has vanishing intensity. The standard propagation 
law for field discontinuities (propagation at the speed of light) is recovered in this limiting 
case. 



It is straightforward to verify that the polar plots of Equation (4.13) are genuine ellipses 
with eccentricity: 



1 + P±A2 ■ 



(4.16) 



Similar observations were made by Boillat, who was however interested in the Born-Infeld 
rather than the Euler-Heisenberg Lagrangian. Note Equation (2.39) in his 1970 paper 



Figure 4.3 visualizes how the effective light cones of (4.10) embed in the background 



geometry. The ^/-dimension is suppressed. Again, four representative cases (the same cases 



used in Figure 4.2) are presented simultaneously for comparison purposes. We note that de 



Oliveira Costa and Perez Bergliaffa [25] have classified effective light cones according to the 
Segre type of the stress-energy tensor for the field. 

Due to birefringence, there are actually two different effective light cones for a given field 



configuration. We did not try to depict both of them in the cases shown in Figures 4.2 and 



4.3 in order to avoid unnecessary clutter. Note however that the difference between the light 
cones of the two polarization states becomes more pronounced at higher intensities. 

As we see from these calculations, the effective light cones are tilted in the direction 
given by the Poynting vector of the field. Based on this observation, we make the following 
conjecture: 

Conjecture 2. For any given carrier field, the corresponding effective light cones tilt into 
the direction of the Poynting vector of the field. Moreover, the higher the field intensity, the 
more pronounced the tilt. 

This conjecture suggests that an optical black hole can form if one contrives to create 
a field, with an inwardly-directed Poynting vector, intense enough to tilt the effective light 
cones all the way to form a trapped surface. 
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X 




-1.0 -0.5 0.0 0.5 1.0 



Figure 4.2: Shock wave fronts in circularly polarized plane wave fields. Four separate 
cases are shown simultaneously for comparison purposes. In a monochromatic circularly 
polarized plane wave of amplitude A propagating in the +z direction, a shock wave initiated 
at the origin is allowed to propagate for a unit t-time. The resulting shock fronts are plotted 



according to Equation (4-13) for four representative cases: (1) P±A = (solid red), (2) 
< P±A < 1 (dashed orange, plotted using P±A = 1/3), (3) P±A = 1 (dot-dashed green), 
and (4) P±A > 1 (dotted blue, plotted using P±A = 3). 

4.5 Distortion of clock readings 

Consider an observer at rest at the origin in t, x, y, z coordinates. Surround the observer 
with clocks, so that in the coordinates these clocks form a sphere S of unit radius. Let 
these clocks be set in such a way that if light travels along null geodesies in the background 
geometry, then the clocks appear to the observer as if precisely synchronized. 

Now assume that the observer is immersed in a plane wave of amplitude A such that 



the effective geometries given by Equation (4.10) pertain. We stipulate that the observer 



sees objects only by way of small disturbances in the plane wave field; the "photons" seen 



by the observer follow null geodesies in the effective geometries (4.10). The readings on 



the stationary clocks at S will no longer appear to be synchronized since the effective null 
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Figure 4.3: Effective light cones in circularly polarized plane wave fields. The transparent 
red, orange, green, and blue cones (outermost to innermost) correspond to the cases plotted 
in Figure In fact. Figure is just the cross section through the plane t = 1. The 
transparent red cone (outermost) corresponds to a standard light cone in the background 
Minkowksi spacetime. The transparent orange, green, and blue cones show how the effective 
light cones embed in the background for plane waves of increasingly intense amplitude. 

geodesies propagate anisotropically witli respect to tlie t, x, y, z coordinates. Moreover, due 
to birefringence, two clock readings may be seen at once. An additional consequence of such 
birefringence effects would be that moving bodies could appear to have double images. 

Since the field disturbances that travel in the direction ip = travel at the usual speed 
of light, the apparent reading of a clock on S as viewed in the direction (f = n will not be 
affected by the effective geometry. By contrast, the other clock readings will be affected. 
One can show that the difference in readings At between a clock viewed at angle if, and 
the unaffected clock at = tt, is given by the formula: 



1 — n(-K 
1 



1 + P^A^ cos^ ip 



P±A^ cos if — 



sin^ ip 



(4.17) 
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Since there are actually two distinct values of P± corresponding to birefringence, there are 
double images. If one sees both polarization states, then two clock readings can be seen. 
In the critical case P±A'^ = 1, the clock reading viewed through ip = is infinitely delayed 
(and not visible), li P±A^ > 1, then the only visible clocks are at angles (f such that 
csc^ > P±A^ . In a special range of cases where P^A^ < 1 < P+A"^, the + polarization 
modes cannot be seen at all when viewed through (y9-angles such that csc^y? < P±A^. 

Since the wave fronts of field disturbances are ellipsoidal in the t, x, z coordinates, one 
might consider reconfiguring S into an ellipsoidal arrangement so that the clocks will appear 
to be synchronized to the observer (provided that the field intensity is kept small enough 
that all points on S are visible to the observer). However, due to birefringence, one would 
only be able to manage the appearance of the clocks as viewed through one polarization 
mode at a time. 



4.6 Hints of an optical black hole? 

In the present chapter, we have found that effective light cones in plane wave fields are 



tilted towards the direction of the Poynting vector of the plane wave (Figure 4.3). Field 
disturbances propagating in the direction of the Poynting travel at the usual speed of light, 
but in other directions there is a drag effect. This drag effect is most pronounced for 
disturbances that propagate in the direction exactly opposite to the Poynting vector. The 
speed of these field disturbances, as measured with respect to the background coordinates, 
is: 

where A is the intensity of the plane wave. 

Comparing effective light cones for plane wave fields of higher and higher intensities as 



in Figure |4.3[ one finds that the light cones become progressively more tilted. A similar 
phenomenon occurs in the geometry of gravitational black holes, where light cones become 
progressively more and more tilted as one approaches the event horizon. At the event 
horizon, the light cones are so tilted that information cannot flow from the event horizon 
to the outside world. We suggest the notion that an optical black hole would form if one 
could increase the intensity of an electromagnetic wave by a sufficiently large amount in a 
localized region of space. 

Though such a field will no longer correspond to a true plane wave, we propose the 
following Gedankenexperiment. Consider an electromagnetic wave that is focusing to a 
point. Let us consider a spherical point-like implosion in which the intensity of the wave 
is assumed to follow the inverse square law. Assuming that the wave front is locally like a 
plane wave, field disturbances that propagate radially outwards would travel at a coordinate 
speed: 

dr r^-P±A^ 
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Here, spherical coordinates (t, r, 9, ip) are implied. Equation (4.19) is calculated by replacing 
A with A/r"^ in Equation (4.18) - in order to take the inverse square law into account. 



Equation (4.19) suggests that a spherical event horizon will form at a critical radius 
P^/^A^^^. That is, within the critical radius, "outgoing" disturbances are not able to escape 
to infinity. 

However, it is not possible to have a nontrivial spherically symmetric electromagnetic 
wave. The fundamental reason for this is that the polarization vectors due to such a field 
configuration would introduce a continuous nowhere-vanishing vector field tangent to the 
2-sphere, thereby contradicting the well-known fact that the 2-sphere is not parallelizable. 

For this reason, we turn our attention to other configurations. Since the cylinder S*^ x M 
is parallelizable, the case of cylindrical collapse can be considered. The next chapter will 
look into this. 

As a tentative calculation for the cylindrical case, replacing A with A/r into Equation 



(4.18) - in order to take the inverse distance law for cylindrical radiation into account - we 



have: 



dr - P±A^ 
di ^ r^ + P±A^'' 



(4.20) 



with cylindrical coordinates {t, r, 6, z) implied, 
horizon at r 



Equation (|4.20|) suggests an effective event 

yield 



1/2 

Pj_ A. Slightly more refined approximations, done in Chapter 



an effective horizon at a radius that is proportional to the square of the intensity A and 
inversely proportional to the frequency. 



Using Equation (4.20) as an estimate for the coordinate velocity of radial null effective 
geodesies which are "outgoing" in a cylindrically imploding wave field, and assuming that 
the ingoing null geodesies propagate at the speed of light, we can use Mathematica to draw 
a graph of the coordinate velocities (see Figure 4.4). More detailed calculations, done in 
Chapter [5| will confirm that this guess is qualitatively on the right track. In Figure -LA, we 



are guessing ingoing rays will propagate at the usual speed of light. This guess is motivated 
by our experience with plane waves; we have seen that field disturbances that propagate 
along with the fiow of a plane wave simply propagate at the usual speed of light. 
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Figure 4.4: Thought experiment. In this graph, we have plotted (in red) our tentative guess 
on how the coordinate velocity v of a radially outgoing light ray will vary with distance r 
in a cylindrically symmetric imploding field. This plot of Equation (4-20) uses P±A^ = 1. 



The horizontal blue line (at v = —1) corresponds to our tentative guess that the ingoing ray 
propagate at the usual speed of light. 
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Chapter 5 
Cylindrical fields 



In Section 4^ we suggested, by means of a crude Gedankenexperiment, that optical black 
holes might arise in the effective geometry of an imploding electromagnetic wave. As we 
have noted, a spherically symmetric implosion is not possible because the 2-sphere S"^ is not 
parallelizable. On the other hand, the cylinder x M by contrast is parallelizable. For 
this reason, in the present chapter we move our focus to cylindrically symmetric fields. The 



main formulas derived in Sections |5.1| and |5.2| have been checked with Mathematical as we 
detail in Appendix [Xj 

Naturally then, we will be working with cylindrical coordinates {t,r,6,z) in which the 
background (Minkowskian) metric g^^, is given by: 



ds' = dt" - dr' - r'de' - dz'. 



The nonvanishing Christoffel symbols are: 

r; 



Or- 



(5.1) 

(5.2) 
(5.3) 



We will define a cylindrically symmetric electromagnetic field as being an A-field whose 
components are functions of the coordinates t and r only: 



At 

A., 



At{t,T), 

Ar{t,r), 
Ae{t,r), 
A,{t,r). 



(5.4) 
(5.5) 
(5.6) 
(5.7) 



Although this is not a necessary condition for making the physical field F^y = d^A^ — (9,^v4^ 
cylindrically symmetric, it is a sufficient one. 



The Euler-Heisenberg field equation (2.21), is: 

„2 



a 
45 



(5.8) 
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This amounts to a system of nonlinear PDEs (as shown exphcitly in Section 5.1 below). 
There are no known general methods for finding exact solutions to such systems. We note 
however, that to first order in a, Equation (5.8) reproduces the Maxwell vacuum equations: 
V^F^" = 0. This means that the familiar exact solutions from Maxwell's theory (whether 
they are cylindrically symmetric or not) are approximations of solutions to Equation (5.8), 
up to first-order in a. 

In Section 5^ we will treat ingoing cylindrical wave solutions from Maxwell's theory 
as approximate first-order solutions to the nonlinear theory. We will plug the Maxwellian 
field into the equations for the effective geometry from Chapter |3] and we will study the 
resulting geometry to second-order in a. As expected, we find that the effective geometry 
is analogous to a black hole. This is evidence for our main conjecture on the existence of 
black hole soliton solutions as we explained in the Introduction. 

Since the Maxwellian solution is only valid up to first-order, one might worry about 
whether it is meaningful to do second-order calculations with it. However, as we show in 
Section 5^, the geometric quantities that we calculate to second-order only depend on the 
first-order (Maxwellian) part of the exact solution, so the approximation is justified. 

In Section 5.3[ we will derive an exact static solution to Equation (5.8). This solution 
corresponds to a constant electric field running in the z-direction together with a magnetic 
field circulating around the 2;-axis. In the linear case, such a field configuration corresponds 
to that of a constant current through an infinitely long straight wire together with a constant 
electric field. With the fields arranged so as to give an inwardly-directed Poynting vector, 
the resulting effective geometry is analogous to that of a black hole. 



5.1 Cylindrical fields of a particular type 

Here we present formulas for the effective metrics and field equations for cylindrically sym- 
metric fields of a particular type. These will be needed in later sections. 

5.1.1 Field tensors 

The particular type of cylindrical field considered here is one in which the t and r-components 
of the A-field vanish. This would be the case, for example, with an elliptically polarized 
imploding cylindrical radiation field in the radiation gauge. 
Let us write: 

At = 0, (5.9) 
Ar = 0, (5.10) 
Ae = u{t,r), (5.11) 
A, = v{t,r). (5.12) 
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The corresponding electromagnetic field tensor F^,^ = dfj,A,y — d^A 



( Fu Ft 



fii/ 



tr 



F'rt Fj-j' Fj.0 F,, 



Fte Ftz \ 

rz 



Fet Fgr Foe Fg^ 
\ Fzt Fzr FzB Fzz / 

/ Q dtu dtv \ 

drU drV 

—dtu ~drU 
y —dtv —drV 



Raising the first index, one gets that = g^^Fx,^ is: 



/ 



^dtu 
\ dtv 



dtu 

—drU 

^dr-U 

'drV 



/ 



dtv \ 

-dr-V 






Raising the second index, one gets that F^^^ = g^'^F'^. is: 



/ 


^dtu —^drU 

\ dtV —drV 



^drU drV 
^ 

/ 



Using (5.13) and (5.15), one gets that the invariant F = F^^F^'^ is: 



The dual tensor F*^^ 



pal3 
( 



is: 






rdrV 

\ -IdrU 



—rdrV -drU \ 

—rdtv \dtu 
rdtv 
-^dtu 



Using (|5.15|) and (|5.17|), one gets that G = F*^F^''' is 



G 



[{dtu){drv) - {dru){dtv)] . 



Raising the indices of the dual, one gets that F*^'^ = g'^^g^'^F^o is: 








r ' 
\ IdrU 






-dtv 



--dtv -dtu 
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5.1.2 Effective metric 



As we saw in Section 3.3, for the Euler-Heisenberg theory, the effective metric g^ji_y is given 
by the inverse of: 



(5.20) 



where: 



A 



224a^ 



± 



495 + 12Fa2 ^ V18225 - 18360Fa2 + 4624^2^4 + 3136G^a^' 



(5.21) 



Henceforth let us write A := A±, keeping in mind that there is a choice of ± involved in the 
calculation of A (this A has nothing to do with the cosmological constant). As we explained 
in Section |3]2, this choice of ± depends on the polarization state of the field disturbance. 
Using Equations (5.14) and (5.15), we get that the only nonvanishing components of the 



substress tensor F^^F^^ are: 



pt^pxt _ 




(5.22) 


pr^pXt ^ 


~^[{dtuM.u)+r\dtv){drv)] 


(5.23) 


pr^pXr ^ 




(5.24) 


pd^pxe ^ 




(5.25) 


pz^pxe ^ 


l^[{drU){drV)-{dtUmv)] 


(5.26) 


pz^pXz _ 




(5.27) 



Plugging our result for the substress tensor into Equation (5.20), we get the effective 
cometric g'^^. Taking the inverse of g'^", we find that the only nonvanishing components of 
the effective metric are (up to a conformal factor k): 



Kgez 



■ ngrt 

Kg„. 

i^gee 
i^gez 

i^gzz 



1 - 



k{drU) 



A{dtu){dru) 



-A{dtV){drV) 

^ Ar'^{drvf - Ar'^{dtvf 
A{dtu){dtv) - A{dru){drv) 
A{druf A{d,uf 

-'-"To o ? 



(5.28) 
(5.29) 

(5.30) 

(5.31) 
(5.32) 

(5.33) 
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Since only the null geodesies are important, we can drop the conformal factor k. Note that 
the effective metric becomes conformally equivalent to the background metric in the limit 
where A — )■ 0. 



5.1.3 Radial null geodesies 



Since the components of the metric tensor only depend on the coordinates t and r, it 
follows that radial null curves (that is, null curves with constant 6 and z coordinates), are 
automatically geodesies. For such a curve we can write: 



= 9u + 2g,r^^+g„i- 



(5.34) 



Solving (5.34) for dr/dt gives: 



dr 
di 



-9t 



gttgr 



{dtu){dru) + r'^{dtv){drv) 



{dtuy + r^{dtvy + 'i 

J[{dtU){drU) + r^{dtV){drV)f + + r-^idtvf + {dtUf) - r^{drVY - {drUY) 

T- 



(5.35) 



We define outgoing geodesies as corresponding to choosing + in (the second line of) Equation 



(5.35) and ingoing geodesies as corresponding to choosing — . 



For radial geodesies of the ingoing type, one gets that to second-order in a: 



dr 
'di 

For the outgoing type: 



dr 
dt 



(5.36) 



out 



(11 ±3)^2 
45^2 



{dtu + druf + r''{dtv + drvf^ + 0{a^). (5.37) 



In Equations (5.36) and (5.37), the choice of ± has to do with the polarization state of the 



disturbance (there is birefringence). With the birefringence averaged out, we have: 

11^2 



and: 



/ dr 
\1t 



dr 
di 



-1 + ^ {{dtu - druy + r^dtv - drvf'^ + O(a^), (5.38) 



out 



IW 

45^2 



(dtu + druf + r^idtv + drvf) +Oia^). 



(5.39) 
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5.1.4 Field equations 

Our next task is to express the field equations as a system of nonlinear PDEs. To this end, 



using (5.15), and introducing u{t,r) := u{t,r)/r, we get that the left-hand side of the field 
equation (5.8) is: 









:^dr{rdru) - \d'^u - 
-dr{rdrv) — d^v 



for u 
for u 
for u 
for u 



t 
r 

e 

z. 



(5.40) 



(5.8) is: 



Using (5.15), (5.16), (5.19) and (5.18), we get that the right-hand side of the field equation 



Here, 



and: 



V 



a 
45 



(4V^(FF^^) + rV^iGF*"")) 






45ri 



for u 
for u 
U for V 



§-^V forz/ 



t 
r 
9 
z. 



(5.41) 



2ru^ [Srd'^u — 3drU — rd^u) — 6u^ + 

r% [6{druY - 2{drvy - h{dtvf - 2{dtu) (dtu + Ardtdru) + 
3r{dtv){dtdrv) + Ar{dru) {3d'^u - d'^u) + r{drv) {Ad^v - 7d^v)] + 

{6{druf - (drv) [7{dtv) {dtu + 2rdtdru) - 3r{dtu){dtdrv)] + 
r{drvf {2dlu + hd'^tu) + {druf {Grd^u - 2rdlu) + 
[dru] [2{drvf - &{dtuf - 8r{dtu){dtdru) + 
{dtv){5dtv + 3rdtdrv) + r{drv) (49> - 7d'^v)] + 
r [(d^^u) {5{dtvf - 2{dtuf) - 7{d',v){dtumv) + 

2{dfu) {3{dtuf + [dtvf) + A{dtu){dtv){d''^v)\ ] , (5.42) 



{2rdlv - 2drV + brd^v) + 
ru [{dtv) (lOSj-u + 3rdtdru) — lAr{dtu){dtdrv) + 
r{drv) {Ad^u - 7d^u) + 2{dru) {2drV + 2rdlv + hrdiv)] + 

{2{drvf - {drv) [2{dtuf - 3r{dtu){dtdru)+ 2{dtv) {dtV + Ardtdrv)] + 
(dru) [3r{dtv){dtdru)- 2{dtu) {2dtV + 7rdtdrv) + 
r{drv) {Adlu - 7d1u)] + {drvf (6r9> - 2rdfv) + 
{druf {QdrV + 2rdlv + brd'^tv) + r [((9» {b{dtuf - 2{dtvf) + 
{dtu){d,v) [Adlu - 7dlu) + 2{dlv) {{dtuf + 3{dtvf)] ] . (5.43) 
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The field equations for a field of the kind specified by Equations (5.9) - (5.12) can thus 
be expressed nonlinear system of PDEs: 



r ' 



[rdj-u) 
rdr-v) 



4oj 
45r' 



u 



4af 
45r^ 



V 



(5.44) 



5.2 Maxwellian approximation 



To first-order in a, in which Maxwell's theory is recovered, the field equations (5.44) become: 



Seeking solutions to (5.45 ) for a monochromatic field of constant frequency w > 0, we express 
the components of the A-field (5.9) - ( 5.12[ ) in the form: 

= Re [S^ exp(iwt)] , (5.46) 



where 5*^ is a function of r only. Note that the ansatz (5.46) implies that d^A^ = — w A^, 
and so (5.45) becomes: 



-dJrdrv) + uj'^v = 0. 



(5.47) 



These equations are of the form: 

1 



n 



(5.48) 



where n is either or 1. Equation (5.48) is a Bessel-type differential equation, having 
solutions of the form (see e.g.. Bowman [30] p. 116): 

ip = ciJn{ujr) + C2Yn{ujr), (5.49) 

where ci and C2 are complex constants, and J„ and Yn denote the nth order Bessel functions 
of the first and second kinds respectively. Monochromatic solutions to (5.45) are therefore 
given by: 



u 

V 



Re [{cuiJi{ujr) + c^s^ilwr)) exp(za;t)] 
Re [(ct,i Jo(wr) + c^2'i^o(w'^)) exp(iu;t)] . 



(5.50) 
(5.51) 



where the Cij are complex constants. 

In order to choose the constants Cij so that one gets radially propagating solutions, 
consider the fact that the graphs of J„ and Yn look like dampened sine and cosine graphs. 
In this sense, combinations such as J„(c<;r) ±iy„(a;r) are like dampened versions of exp(ici;r). 
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So, after taking the real part, the functions (J„(a;r) ± iYn{ur)) exp{iojt) describe radially 
propagating waves (which can be either ingoing or outgoing depending on the choice of ±). 
Accordingly, for an elliptically polarized ingoing cylindrical wave, we write: 



u(t, r) 



r ■ Re 



U 



{Ji{ur) + iYi{ur)) exp{iut) 



and: 



v(t, r) 



Ur 

CO 



Re 
V 

OJ 



Jiiojr) coiiiojt) — Yi{ojr) sm{ojt) 



(5.52) 



V 



OJ 



{jQ{ojr) + iYQ{ojr)) exp(za;t) 



Jo(u;r) cos((X't) — Y^iur) sm{ojt) 



(5.53) 



where U and V are real-valued constants (not to be confused with the functions U and V 
as defined by ( |5.42| ) and ( |5.43[ )). 

Our task is to study the effective metric corresponding to this wave field. More specif- 
ically, we want to have a look at the effective radial null geodesies using (5.38) and (5.39) 
to calculate dr/dt to second-order in a, where we use the Maxwellian solution to evaluate 
the field variables. We claim that, up to second-order in a, the calculation of dr/dt only 
depends on the first-order (Maxwellian) part of the exact solution to (5.44). In other words: 

Theorem 8. Maxwellian approximations for dr/dt (along radial null geodesies in the effec- 
tive geometry) are accurate up to second-order. 

Proof. Suppose that we have an exact solution {u,v) = {se/r,Sz) for the nonlinear system 
(5.44). In the limit — > 0, this solution becomes a solution to (5.45). So expanding the 
exact solution as a power series in a would yield sg = mg + O(a^) and = rriz + O(a^), 
where («,!>) = {mg/r,mz) is an exact solution to (5.45). (Note: a priori the series may only 
be asymptotic.) Using Equation (5.39): 



/dr\ 
\dt/ 



1 - 
1 
1 - 



45r2 
45r2 
45r2 



{dtSe + drsgf + r\dts, + drSzf) + 0{a^) 
{dtmg + drmg + 0{a^)f + r'^{dtm, + drm, + 0{a^)f 
{dtmg + dringf + r'^{dtmz + <9,.m^))2j + 0{a^). 



Similar calculations can be done using Equations (5.36) - (5.38). 



O(a^) 
(5.54) 

□ 



Proceeding now, by specializing Equations (5.38) and (5.39) to the Maxwellian solution 
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(5.52) and (5.53), we find that the radial null geodesies are described by: 



dr 
dt 



-1 



45 



{Jo{ujr) + Yi{ujr)) cos{ujt) + {Ji{uir) — 10(1^^)) sin(u;t) 



(10(1^^) — Jii^r)) cos{ujt) + {Yi{ijjr) + Jq^ut)) sin(a;t) 



+ 



(5.55) 



and: 

/dr 
\dt 



1 - 



1^ 



(Yi^ur) — Jo(u;r)) cos{ut) + {Jiiuor) + 1^0(1^^)) sin(a;t) 



1 2 



+ 



(Ji(u;r) + YQ{ur)) cos{ujt) + (Jo(wr) - Yi(wr)) sin(wt) \ + 0{a 



(5.56) 



We note that the oscillatory terms involving trigonometric functions of t disappear in the case 
of circular polarization (where U = V). One might have expected this out of consideration 
of the fact that, as we saw in Chapter |4| a similar simplification occurs in the effective 
geometry of circularly polarized plane waves. In fact, the stress-energy tensor (at least, as 
computed using the Maxwellian Lagrangian L = —F/A) for the field given by Equations 
(5.52) and (5.53) does not have any oscillatory terms in the case where U = V. 

Henceforth, let us assume that the wave is circularly polarized, with U = V =: A. In 
this case, we get that the effective radial geodesies are described by: 



/dr 
\di 



-1 + 



45 



■Kur 



+ Joicorf + Ji{ujrf + Y^{urf + Yi{urf + 0{a^) 



(5.57) 



and: 



dr 
di 



1 - 



out 



45 



Tcur 



+ Uurf + J,{urf + Y^iurf + Y^{urf + O(a^) 



(5.58) 



For simplicity, we are using the formulas in which the birefringence is averaged out. To 
recover the birefringence, replace the factor lla^ with (11 ± 3)q;^. 

We note that, for large x, one has the approximations (Arfken and Weber [21] p. 718): 



JJx) 




cos 



X — \ n + 



and: 



YJx) 



sm 



TtX 



X 



n + 



(5.59) 



(5.60) 
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Consequently, for large x\ 



Mxf + Ji{xf + Y^{xf + Yi{xf 



cos 



X 



H cos 



H sm 



H sm 



X — 



1 + 



1 
1 



X - 1 + 



(5.61) 



In fact, one could have guessed at Equation (5.61 ) using the following idea. Far from the 



origin, the cylindrical wave should look like a plane wave and we know that field disturbances 
in a plane wave, when they travel along the same direction as the plane wave (the Poynting 

vector), travel at the usual speed of light. Hence one should have 



—1 in the limit 



where ur is large, and this implies (5.61). 



So in the limit where the quantity lot is large, we have, to second-order in a: 

dr 
dt 



and: 



/dr 
\dt 



out 



Equation (5.63) suggests that within radii r < Vc, where: 



4:57100 



(5.62) 



(5.63) 



(5.64) 



even the "outward" geodesies are compelled to fall inward. Hence the critical radius rc is 
the event horizon of a black hole. (Note: we have only checked this for outward geodesies 
in the radial direction.) 



Since (5.64) was derived assuming that ur is large, we only expect this approximation to 
hold in the limit of very large A^. (We could also mention that, due to birefringence, there 
are actually two event horizons. If we had taken this into account in the above, then we 
would have estimated the critical radii as occurring at ^ (88 ± 2A)a^A'^/{A5nu), where 
the ± depends on the polarization of the disturbance.) 



Plotting Equations (5.62) and (5.63) on the same graph (Figure 5.1), we can compare to 
our earlier naive guess of (4.20). We note that there are substantial quantitative differences 



between our initial guess and our slightly more refined calculation, but the qualitative picture 
is basically the same. 
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Figure 5.1: In this graph, the coordinate velocities of effective null geodesies are plotted 



using the asymptotic approximations (5.62) and (5.63), which assume that the quantity lot 
is large. The horizontal blue line at dr/dt = —1 corresponds to the ingoing geodesies, and 



the red curve corresponds to the "outgoing" radial geodesies. Compare to Figure 4-4 



Treating tlie asymptotic approximations (5.62) and (5.63) as ordinary differential equa- 



tions, and solving them by integration, we obtain approximate equations for the radial null 
geodesies. Specifically: 



t 



-r + ro. 



(5.65) 



for the ingoing geodesies (ro := the radial coordinate of the geodesic when t = 0), and: 



r + Tc ln(r — Vc) — tq — Vc ln(ro — Vc) if ro > Tc 
r + Tc ln(rc — r) — Vq — r^. ln(rc — ro) if ro < r^ 



(5.66) 



for the "outgoing" geodesies. An outgoing radial null geodesic that initiates from r 



would just remain there. Using Mathematical we have plotted Equations (5.65) and (5.66) 



for a few values of ro. The resulting plot is shown in Figure |5.2 Note that the effective 
light cones are tilted in towards the origin, just like in the situation of gravitational black 
holes. 

In fact, using the numerical integration capabilities of Mathematica, we can make space- 



time diagrams for the effective null geodesies, as described by Equations (5.57) and (5.58) 
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1.0 



0.5 - 



event horizon a) r = r, 



outgoing geodesic with rg > 



-0.5 



-1.0 - 



-1.5 - 



-2.0 



0.5 



1.5 



2.0 



2.5 



3.0 



ingoing~^odesic = 



ingoing-geodesic =0 



Figure 5.2: In this graph, effective null geodesies are plotted using Equations (5.65) and 



(5.66). Equations (5.65) and (5.66) are themselves based on the asymptotic approximations 



(5.62) and (5. 65), which assumes that the quantity oor is large. 



to second-order in a, without recourse to the asymptotic approximations (5.62) and (5.63). 



These diagrams, shown in Figures 5.3 - 5.5, are similar to Figures 5.1 - 5.2 Again we find 



that the effective geometry contains a black hole. That is, the effective light cones are tilted 
towards the origin, and there is an effective event horizon. In making the plots for Figures 

1 and we have chosen such that SSa'^A'^/ (457ru;) 



5.3 - 5.5 



we have set a; = 1 and we have chosen A'' such that 88a^y4^/(457ru;) = 1. 
As we see from these plots, the radially outgoing rays are significantly slowed down near 
the critical radius. This means that if we slowly move a clock radially inwards, then an 
observer at infinity would see it ticking at a progressively slower rate. Actually, due to 
birefringence, the situation is even more complicated since there will also be double images, 
but we are ignoring birefringent effects for now. When the clock reaches the critical radius, 
its light rays will not travel beyond the critical radius, and the clock will no longer be visible 
from the outside. 



We note that, according to the second-order approximations (5.57) and (5.58), there is 



a small radius within the event horizon where the ingoing geodesies are brought to zero 
coordinate velocity. Thereby the ingoing geodesies coming in from infinity do not penetrate 
all the way to the origin, but instead are blocked by an "inner horizon" at r = r/ (see 
Figure (5.4)). Also, we note that null geodesies exceed the usual speed of light, as viewed 
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Figure 5.3: In this graph, the coordinate velocities of effective null geodesies are plotted 
using Equations (5.51) and (5.58) up to to second-order in a. In plotting this graph, we 
have set u = 1 and we have chosen A? so that the quantity 88a;^A^/(457ru;) (our crude 
estimate for the effective horizon radius) is unity. The blue curve corresponds to the ingoing 
geodesies, and the red curve corresponds to the "outgoing" radial geodesies. Compare to 
Figure \5.1 



in tlie background coordinates. Some of tlie phenomena sliown in Figure 5.4, especially at 
the smaller radii, may be mere artifacts of the approximation. We note however that the 
superluminal photons, if such exist, will not violate causality if the effective spacetime which 
they propagate is a causal spacetime. 

Note also that there is a radius between the inner and outer horizons where the "ingoing" 
and "outgoing" geodesies cannot be locally distinguished. At this special radius, the ingoing 
and outgoing radial geodesies travel in the same direction at the same velocity, so the 
effective light cone is degenerate at this radius. 

We remark that if one were to take u as negative, which amounts to turning our ingoing 
wave into an outgoing cylindrical wave, then one would obtain an effective spacetime which 
contains an optical white hole rather than a black hole. 
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Figure 5.4: In this graph, effective null geodesies are plotted by numerically integrating 
Equations (5.51) and (5.58). We have set u = 1 and we have chosen A? so that the 



quantity SSa"^ / (Abnu) is unity. Compare to Figure 5.2 
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5.3 A static exact solution 



In the case of a field with 2;-polarization, u = 0, the field equations (5.44) reduce to a single 
nonlinear PDE: 



drV 



r 



to" 
45r 



[{dj-vY — {dtv){drv){dtv + Ardtdrv) 



+r{dtvf{?>d'^v - dlv) + r{drvf{?,dlv - 5^^)] . (5.67) 



A particularly interesting case in which Equation (5.67) can be solved exactly is that of a 
field where dtv = E = constant. In this case (5.67) reads: 



d'^v 1 dv 
dr'^ r dr 



8a' 
45r 



/ dv 
\ dr 



/dv^^ 
[d^J 



rE^ 



dr'^ 



+ 3r 



dv\ 
dr J 



dr'^ 



Introducing the function B(r) := dv/dr, Equation (5.68) becomes: 



dr r 



45r 



E^B 



fdB\ 



rE^\ — \+ 3rB' 



fdB\ 



This can be rearranged into: 

dB 
dr 



Sa'^B^ - Sa^E'^B - 455 



(5.68) 



(5.69) 



(5.70) 



j(24a252 _ 3c^2^2 _ 45) 

Equation (5.70) is solvable by integration. The general solution is given through the relation: 

8«2 



B + (E^B - B^) 

45 ^ ^ 



k 



(5.71) 



where k is an arbitrary constant. A sketch of the graph of Equation (5.71) is shown in 
Figure |5.5[ There are three asymptotic values of i? as r — )■ oo, namely: B = 0, and 



Let 



B = ±^E^ + A5/{8a'^] 

Note that dB / dr has a singularity when r = and when B = ± a/ E'^/3 + 15/(8a;2 
us define Bg := ^/Ey3~+T5Ji8oFj , and let denote the radius where B'^ = B^. 

We remark that (5.71 ) corresponds physically to a constant electric field E directed along 
the ±z-direction (± depending on where E is positive or negative, respectively), together 
with a magnetic field B, which in the Maxwellian limit — )■ 0, would be produced by a 
constant current 2iTk running along the ±2;-direction (± depending on whether k is positive 
or negative, respectively). The magnetic field circulates around counterclockwise around 
the the z-axis if B is positive, clockwise if negative. In other words. Equation (5.71) refines 
the familiar undergraduate physics formula B 



k/r. 

Using Equations (5.16) and (5.18), we get that the F and G invariants are: 

F = 2{B^-E'), 



(5.72) 
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B 



Figure 5.5: A sketch of the graph of Equation (5.71) for k > 0. 



and: 



G 



(5.73) 



So Equation (3.44) gives: 



1 

A 



495 + 24^2(52 - E^) T 135 - 136a^{B 



E^ 



224^2 



(5.74) 



where the choice of ± depends on the polarization state of the field disturbance. We shall 
call the polarization corresponding to choosing + in Equation (5.74) the "(+) polarization 



state," and we call the other state the "(— ) polarization state." (Note that our naming 
schemes for the polarization states are always ad hoc and the naming scheme in the preset 



section is not meant to be consistent with e.g., Section 4.1 or Appendix A 



Using Equations (5.28) - (5.33), we get that the only nonzero components of the effective 



metric are, up to a conformal factor: 

9tt 

Qtr = 9rt 



l-AB^ 

-AEB 

-1-AE^ 

-r^ + Ar^ [B' - E^) 
-1. 



(5.75) 
(5.76) 
(5.77) 
(5.78) 
(5.79) 
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Using Equation (5.35), we get that the radial null geodesies are given by: 

-EBt 



dr 
~dt 



A2 A 



\ (B' - E^) 



(5.80) 



The outgoing radial geodesies correspond to choosing + in Equation (5.80), and the ingoing 
geodesies correspond to choosing — . 



5.3.1 Proof of the main theorem 



Since Equation (5.71) gives a cubic equation in B, there are three branches giving i? as a 



real function of r (see the Table below). 



branch 


range of 


I 


to 52 


II 


Bs to 


III 


to 00 



Table 5.1: The three branches of B classified according to their ranges. 
Let us consider the case where the field is prescribed by branch I, with k > 0. This field 



is defined only for r > r^. Note that Equation (5.70) can be rewritten as: 

?2N 



dB 
dr 



B (352 - B' 



Bl) 



(5.81) 



3r (fi2 

Thus we have that i? is a strictly decreasing function of r. As shown below, both effective 
geometries for this field contain black holes, if \J~^^ < E < \J~^^- This section constitutes 
proof of Theorem [T] from Chapter [Tj 

Theorem 9. For E"^ > ^^2, the effective geometry corresponding to the (+) polarization 
state contains a black hole (if E > 0), or a white hole (if E < 0), with the event horizon at 
r = r^. 



Proof. For the (+) polarization state, with E'^ > 45/(34q;2), Equation (5.74) gives: 

45 



1 

A 



16a 



At r = r,: 



1 

A 



Bi 



(5.82) 



(5.83) 



In fact, since 1/A > 1/A|r=j.^, we have that A > for all r > r^. 
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Equation (5.80) yields: 

-EB- 



dr 
di 



out, r>rs 



-EB+J{^-1 {B^ - E-)Y - (if. - i (5^ - E^)) {B^ - E^) 



E^ + 



45 
16a2 



I (52 - E^) 



(5.84) 



For r > Vg, the numerator in Equation (5.84) is positive by virtue of the fact that B^ < B^ 
in the region r > r^. Since A > 0, the denominator is also positive. So dr /dt\ont, r>rs is 
positive in the region r > r<j. 

For ingoing radial null geodesies, we have: 



dr 
dt 



-EB 



45 



I (52 - E^)Y - (if. - - E^)) - E^) 



in, r>rs 



E^ + 



45 



I (52 - E2) 



(5.85) 



and dr/dt\i,^^ r>rs is negative since 52 < 52 in the region r > r^. 
At r = rg, we get: 



B,{\E\-E) 



out, r = rs 



E^ + B 



2 ' 

s 



(5.86) 



and: 



dr 
'dt 



B,{\E\+E) 



(5.87) 



in, r = r. 



If ii^ is positive, then dr/dt 



E^ + 52 ■ 

0; radial outgoing null geodesies at r = are trapped. 

= 0; radial ingoing geodesies cannot 



in, r = Ts 



out, r = r. 

On the other hand, if E is negative, then dr/dt 

reach r = Ts from r > Tg (the outside). 

It remains to be shown that the nonradial curves are trapped at r = r^. To this end, 
suppose that we have an arbitrary null curve in the effective spacetime. We write: 



dv 

= gtt + '^9tr^ +grr [ 



V dt. 



fdd 

Vdt 



+ 9z 



dz\' 
dt / 



(5.88) 



with gi^ii, given by Equations (5.75) - (5.79). Then: 



dr 
dt 



-9tr ± Jgtr - gttgrr - ^rr (^eS (f )^ + )' 



(5.89) 
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We claim that: 



-Qtr + V Qtr - 9tt9rr < ^ <- ^dtr - a/^ 

- dt - 



9tt9r 



9r 



9r 



(5.90) 



That is, an arbitrary null curve cannot climb up to larger radii faster than a radially outward 
geodesic, and cannot fall down to smaller radii faster than a radially inward geodesic. In 
other words, nonradial null curves are trapped if the radial null geodesies are trapped. 

Note that in order to prove the claim, it suffices to show that Qrr < and gge < 0. 

We get that Qrr {= ~^ " Ai?^) is negative since A > 0. 

To get gee = -r^ + Ar^iB^ - E^) < 0, it suffices to show that B^ - E^ < 1/A. To this 
end, note that: 



B^ -E^ 



15 

< 



E' 



2 ■ 



Muliplying both sides of (5.91) by 3/2, we get: 

Ub^-e-)< 



(5.91) 



(5.92) 



so: 



B^ - E' 



< 



45 
16a2 

1 

A' 



[B^ - E') 



(5.93) 

□ 



Theorem 10. For < E^ < the effective geometry for the (— ) polarization state 
contains a black hole (if E > 0), or a white-hole (if E < 0), with an effective event horizon 
at r = Tc such that r > r^. In fact, = 9k\/5/ (7aE'^\/18 — 8a^i?^) . 

polarization state, with E^ > 45/(34a;^), Equation (5.74) gives: 

1 _ 45 
A ~ 



Proof. For the 



(5.94) 



Since E'^ < 9/(4a;^), we have that A > 0. 
At r = r^: 

1 

A 



495 - 80«2^2 
168^2 • 



Moreover, for outgoing radial null geodesies. Equation (5.80) gives: 



dr 
di 



'-168a^E + 4av^30(99 - 16a^E^) 
495 + 88^2^2 



5. 



(5.95) 



(5.96) 
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which is nee; ative a E > +^/A5/{3Aa^. This means that the outgoing radial geodesies, 
initiated from r = r^, are compelled to fall down to smaller radii. At the other extreme 
note that: 



oo 



dr 
dt 



1 



out, r=oo 



> 0. 



(5.97) 



Hence, by the intermediate-value-theorem, there is a radius rc, which is greater than and 



less than oo, where dr/dt 



out, r= 



(black hole event horizon at if E > y^A5j{3Act^) . 



In fact, the critical radius Vc is unique and we can calculate it. If we set the left hand 



side of Equation (5.80) equal to 0, and solve for B, then we find that there is only one 



real-valued positive solution, namely: 



Be 



45 - 20^2^2 



(5.98) 



Thereby, using (5.71), we get: 



5, + ^ (E^B, - Bl 



9kV5 



7aEWl8 - 8a^E^ 



(5.99) 



Note that dr/dt\o^t changes sign at r = Tc since dr / dt\out,r=rs is negative and dr / dt\ out, r=oo 
is positive. That is, any outgoing geodesic in the region r < is compelled to fall inward 
to smaller r; any outgoing geodesic in the region r > r^ will escape to larger r. 



Next we consider the ingoing null geodesies. For these. Equation (5.80) gives 



dr 
~dt 



IQSa^E - 4av^30(99 - 16«2^2)\ 

rrr^ — — -D« 



495 + 88^2^2 



(5.100) 



This is positive li E < -^A^/{Ma'^). At r = cxD, we have: 

1 



dr 
di 



m, r=oo 



Vl + AE2 



< 0. 



Thus, if E is negative, there is a radius r^, between and oo, where dr/dt 



(5.101) 







(white hole event horizon at r^). The quantity dr/dt changes sign at Tc in such a way 

in 

that ingoing geodesies issuing from the region r < rc will be compelled to escape outward 
to larger r, and ingoing geodesies issuing from r > rc will fall inward to smaller r. 

As in the proof of the previous theorem, we get that the nonradial geodesies are trapped 
by the event horizon by showing that gge < and ^Vr < 0. The fact that A > gives Qrr < 0. 
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As before, to get gee < 0, it suffices to show that B — E < 1/ A. To this end, note that: 



< 



45 



(5.102) 



Multiplying both sides of (5.102) by 2/7, we get: 

-(B'-E') <^. 
7 ^ ' 28a2 ■ 



(5.103) 



so: 



28a2 ^ 7 V ) 



1 

A' 



(5.104) 

□ 
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Appendix A 
Mathematica notebook 



The purpose of this appendix is to show how Mathematica (version 8) can be used to check 
or carry out the calculations in 5A_ and |5.2 Our implementation of tensor calculus in 



Mathematica is modeled on applications found elsewhere, such Miiller and Grave [32] and 
the downloadable notebooks of Parker [55] . 

Let us begin by clearing out the variables that will be used: 

Clear[coord, t, r, 9, z, i,j, k, /, u, v, metric,cometric, aflfine, Afield, faraday, faradayl, faraday2, 
Fspecialcase, maxwell, Gspecialcase,maxwelll, maxwell2, substress, ecometric. A, emetric, 
simplifiedemetric, a, lambdaplus, lambdaminus, radAl, radA2, radBl, radB2, CDfaraday2, 
s, o, Ffaraday2, CDFfaraday2, Gmaxwell2, CDGmaxwell2, righthandside, £/, V, a;, A] 

Next, we specify the coordinate system (cylindrical coordinates) and the background 
metric (Minkowski spacetime): 

coord: =coord = {t, r, 6, z} 

metric:=metric = {{1, 0, 0, 0}, {0, -1, 0, 0}, {0, 0, -r^, 0} , {0, 0, 0, -1}} 

The background cometric is computed by entering: 
cometric:=cometric = Inverse [metric] 

The Christoffel symbols (for the background metric) are calculated by entering: 
aflSne:=affine = FullSimplify[ 

Table[ |Siim[(cometric[[A:,/]])* (D[metric[[i,Z]],coord[[j]]]+ 
£)[metric[[7, 1]], coord[[i]]]— £)[metric[[i, j]] , coord[[Z]]]), 
{Z,l,4}],{fc,l,4},{z,l,4},{i,l,4}]] 
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A.l The field tensors 



Our first task is to check the results given in Section 5.1.1 To this end, we input the A-field 



that we wish to study (coinciding with Equations (5.9) - (5.12)): 
Afield:=Afield = {0, 0, u[t, r],v[t, r]} 

Next we computer the field tensor F^^, called "faraday:" 
faraday:= faraday = 

Table[D[Afield[b]], coord[[z]]]- £)[Afield[[z]], coord[b]]], {i, 1, 4}, {j, 1, 4}] 

The components of F^j^ are displayed by entering: 
faraday/ /MatrbcForm 



\ -t,(i.o)[t,r] -f(°'i)[t,r] / 



(A.l) 



In Mathematical r] denotes dtu^ and M'^°'^^[t,r] denotes drU^ etc. This output agrees 



with Equation (5.13) 



Raising the first index, we get -F^V- Call this "faradayl:" 
faraday 1:= faradayl = FullSimplify[ 

Table[Sum[cometric[[«, fc]]faraday[[fc, 1, 4}], {«, 1, 4}, {j, 1, 4}]] 



Displaying the components of F^^y in matrix form, as in Equation (5.14): 
faradayl/ /MatrixForm 



/ 











^(1,0) 
_y{0,l) 



,(1.0) 



[t,r] 



,(0.1) 



[t,r] 



^;(l,0) 
_y(0,l) 



[t,r] 













Raising the second index, we get F^^^ ( "faraday2" ) : 



faraday2:= faraday2 = FullSimplify[ 

Table[Sum[cometric[[A;,j]]faradayl [[«,/!;]],{/;;, 1,4}], {i, 1,4}, {j, 1,4}]] 



(A.2) 
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As in Equation (5.15), we have: 



faraday2/ /MatrixForm 










M(l'0)[f,r] 






M(0'l)[t,r] 



We get that the F-invariant is, in agreement with (5.16): 



Fspecialcase = Simplify[Sum[faraday[[i,j]]faraday2[[«, j]],{i, 1,4}, {j, 1,4}]] 



For the dual tensor F* ("maxwell":) 

maxwell:=maxwell = FiillSimplify[Table [|Sqrt[— Det[metric]] 
Sum[Signature[{i,j, fc, Z}]faraday2[[i, j]],{i, 1, 4}, {j, 1, 4}], {fc, 1, 4}, {/, 1, 4}], i 



Displaying F* matrix, as in Equation (^5.17|: 

maxwell/ /MatrixForm 



/ 

-ri;^"'"^|t,r 

\ [f,r] M(^-°)[t,-r] Q 





/ 



The G-invariant is, in agreement with Equation (5.18): 



Gspecialcase = Sum[maxwell[[i, j]]faraday2[[i, j]],{i, 1,4}, {j, 1,4}] 



4^,(0,1) ^]^(i,o) ^] 4^(0,1) ^]^;{i,o) ^] 



To get F*^,,, enter: 



maxwelll:= maxwelll = FullSimplify[ 
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Table[Sum[cometric[[i, fc]]maxwell[[/i;, j]\,{k, 1,4}], {i, 1, 4}, {j, 1, 4}]] 

For F*'^'^, enter: 
maxwell2:= maxwell2 = FullSimplify[ 

Table[Sum[cometric[[fc, j]]maxwelll[[i, k]],{k, 1,4}], {i, 1,4}, {j, 1,4}]] 



The components of F*^'^, as in (5.19), are: 
maxwell2/ /MatrixForm 



/ 



V - 





t)(0'i)[f,r] 
(0'i)[i,r] 






[t,r] 



i;(0-l)[f,r] 

^(1^0) [t,r] 
r 











(A.7) 



A. 2 The effective metric coefficients 



Now let us check Section 5.1.2 To this end, we calculate the substress tensor F^\F'^'^: 



substress:= substress = FullSimplify[ 

Table[Sum[faradayl[[i, fc]]faraday2[[fc,i]],{/c, 1,4}], {%, 1,4}, {j, 1,4}]] 

The nonzero components of the substress are displayed upon entering the lines: 

Do[If[UnsaineQ[substress[[i, j]], 0], 

CellPrint[ 

Display Form [ 

RowBox[ 

{SubscriptBox[SuperscriptBox["F" , coord[[i]]], 
"A" ] , SuperscriptBoxf "F" , 

RowBox[{"A'',coord[b]]}]],"=", 
substress[[i, j]]}]] 
]],{i,l,4},{j,l,4}] 
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P0 T^Xe _ »(°'i)[f,r]2-«(i.«)[f,r]2 
XJ^ — ,,4 

pe^pXz _ m(0>i) [t,r]ti(°-^) [t,r]-M(i-») [t,r]t)(i'») [f ,r] 



\e _ m(O.i) [t,r]i;(°-^) [t,-r]-M(i-0) [t,r]t)(i'") [f ,r] 



r2 



The output above agrees with Equations (5.22) - (5.27). 
To calculate the effective cometric, enter: 

ecometric:=ecometric = cometric + A substress 



Here, A stands for A± in Equation (3.44). The effective metric is calculated by entering: 

emetric:=emetric = Inverse [ecometric] 

Let us multiply the effective metric by a certain conformal factor; this choice of conformal 
factor considerably simplifies the effective metric coefficients: 

simplifiedemetric:= simplifiedemetric = 

FullSimpUfy[ - (((-1 + A [t,r] 2) (r^ + A u^'^^^^rf) - 

2A2m(o.i) [t , r]?;(0'i) [t , r]^^!'") [t, r] % r] - 

r^A t;(i>")[t,r]2 + A u^'^^^^rf (l + A v^'^^^.rf)) /r2)emetric] 

The nonzero components of the conformally rescaled effective metric are displayed upon 
entering the lines: 

Do[If [UnsameQ [simplifiedemetric[[i, j]] , 0] , 
CellPrint [Display Form[ 
RowBox[{coord[[i]], coord[[7]], "-comp", 
"=" , simplifiedemetric[[i, j]]}]] 
]],{*, 1, 4}, {j, 1,4}] 

tt - comp = 1 - A!i^!^jM! _ A^(o,i)[t^ rf 

tr - comp = - ^"'°'''[^;t'''°''*'^] - At;(o.i)[t,r]t;(i'0)[t,r] 

rt - comp = - ^"'°'''[*;t'''°' - At;(o,i)[t^ r]v'^^fi^[t, r] 

rr — comp = — 1 Af ' [r, rj 

ee - comp = r2 (-1 + Ai;(0'i)[t,r]2 - Av^'^'^^[t, r]^) 

Oz - comp = -Am(°'1) [t, r]t;(°'i) [t, r] + Am(^'°) [t, r]f [t, r] 

^6^ - comp = -Am(°'^) [t, r]w(°'^) [t, r] + Am(^'°) [t, rjf (^'^^ [t, r] 

— comp = 



This output agrees with Equations (5.28) - (5.33) 
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A. 3 The radial null geodesies 



Now we check Equations (5.36) and (5.37) from Section 5.1.3, which concerns radial null 



geodesies in the effective geometry. To this end, we will need to let Mathematica compute 



the values A± from Equation (3.44). The value A+ is designated "lambdaplus," and A_ is 
designated "lambdaminus:" 

lambdaplus: = lambdaplus = (224q;^)/(495 + 12Fspecialcase a^— 

Sqrt [18225 - 18360Fspecialcase 4624Fspecialcase^a^ + 3136Gspecialcase^a*] ) 



lambdaminus: = lambdaminus = (224q;^)/(495 + 12Fspecialcase 

Sqrt [18225 - 18360Fspecialcase aH 4624Fspecialcase^a^ + 3136Gspecialcase^a*] ) 

The choice of ± in the calculation of A± depends on the polarization state of the photon. 
We call these the + and — polarization states (corresponding to A+ and A_, respectively). 



Using Equation (5.35), we get that for outgoing geodesies in the + polarization state. 



dr/dt expanded as a series in a is: 
radAl = PuUSimplify [Series [ 

i-D[u[t, r], t]D[u[t, r], r]- r^D[v[t, r],t]D[v[t, r], r]+ 
Sqrt[ (D[u[t, r],t]D[u[t, r],r]+ r^D[v[t, r], t]D[v[t, r], r]f + 



( 

Vlambd 



lambdaplus 

,2 



+ 
+ 



lambdaplus 

{a,0,3}],r2>0] 



D[u[t,r],tr + r^D[v[t,r],tf)])/ 
D[u[t,r],t]^ + r^D[v[t,r],tr), 



For outgoing geodesies in the — polarization state: 



radA2 = FiillSimplify[Series[ 

r], r], r]- r2Z)[?;[t, r], t]£>[?;[t, r], r]+ 

Sqrt[ {D[u[t, r],t]D[u[t, r], r]+ r^D[v[t, r], t]D[v[t, r], r]f + 
(i^^b=i - D[u[t,r],r]'-r'D[v[t,r],rr) 

^ lambdaminus 

(i^^b^ + r], tr + r'D[v[t, r],tr) , 

{a,0,3}],r2>0] 
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1 - 



45^2 



Outputs (A. 8) and (A.9) imply Equation (5.37) 



For the ingoing geodesies with + polarization: 
radBl = FullSimpUfy [Series [ 

{-D\u[t, r], t]D\u[t, r], r]- r''D[v[t, rlt]D\v[t, r], r]- 
Sqrt[ {D[u[t,r\,t]D[u[t,r\,r\+ r'^D[v{t,r\,t]D[v[t,r\,r\f + 



f lambdaplus 
^ lambdaplus 



lambdaplus 

{a,0,3}],r2>0] 



- 1 + 



1 



45r2 



14 ( [t , r] - [t, r] ) ' + [t, r] - v^^^^^ [t,r])^) a" + O [a] 



(A. 10) 



For the ingoing geodesies with — polarization: 
radB2 = FullSimpUfy[Series[ 

{-D\u\t, r], t]D\u\t, r], r]- r'^Dlvlt, rlt]D\v[t, r], r]- 
Sqrt[ {D[u[t,r\,t]D[u[t,r\,r\-^ r^D[v\t,rlt]D[v[t,r\,r\f + 

(i^b^ - rl rf - T'D[v[t, r], rf) 



lambdaminus 
^2 



+ 
+ 



Y lambdaminus 

{a,0,3}],r2>0] 



£)Mt,r],tp + r2D[t;[t,r],t]2)])/ 
£)Mt,r],t]2 + r2£)[t,[t,r],t]2), 



(A.ll) 



Outputs (A. 10) and (A.ll) imply Equation (5.36) 



A. 4 The field equations 



Our next task is to derive the nonlinear PDEs (5.44) which arise from the Euler-Heisenberg 



field equation (2.21) together with the cylindrical field ansatz presently under consideration. 
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Since we wish to express these PDEs in terms of the functions v and u = u/r, we will 
back and re-enter the field tensor, so that in the present section, u[t,r] should be read 
" u[t,r] ". To this end, we will clear out and recompute the relevant variables: 

Clear[Afield, faraday, faradayl, farada3r2, Fspecialcase, 

maxwell, Gspecialcase, maxwelll,maxwell2] 

We re-enter the 74-field as: 

Afield:=Afield - {0, 0, r u[t, r], v[t, r]} 

Recalculating Fn,^: 

faraday:= faraday = Table[£)[Afield[[7]],coord[[«]]]— 
I>[Afield[[i]],coord[b]]], {i, 1,4}, {j, 1,4}] 

I/" 

faradayl:= faradayl = FullSimplify[Table[Siim[cometric[[i, fc]]faraday[[A;, j]], 
{fe, 1,4}], 0,1,4}, 0,1,4}]] 

faraday2:= faraday2 = FullSiinplify[Table[Siim[cometric[[fc, j]]faradayl[[i, fc]], 
{A;, 1, 4}], 0,1, 4}, {i, 1,4}]] 

F: 

Fspecialcase:= Fspecialcase — Simplify[Siim[faraday[[i, j]]faraday2[[i,_7]], 
{z,l,4},{j,l,4}]] 

fr* ■ 

maxwell:=maxwell = FullSimplify[Table [lSqrt[— Det[metric]] 

Siiiii[Signature[{i, j, fc, Z}]faraday2[Ii, j]],{i, 1, 4}, {j, 1, 4}], {fc, 1, 4}, {/, 1, 4}], r > 0] 

G: 

Gspecialcase: =Gspecialcase — Siiin[maxwell[[i,j]]faraday2[[«,^]], 
{i,l,4},{i,l,4}] 

V 

maxwelll:= maxwelll — FullSimplify[ 
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Table[Sum[cometric[[i, fc]]maxwell[[/i;, j]\,{k, 1,4}], {i, 1, 4}, {j, 1, 4}]] 



maxwell2:= maxwell2 = FullSimplify[ 

Table[Sum[cometric[[fc, j]]maxwelll[[i, k]],{k, 1,4}], {i, 1,4}, {j, 1,4}]] 

Next, we calculate and display the covariant derivative V ^F^'^ by entering: 
CDfaraday2:= CDfaraday2 = FullSimplify[ 

Table[Sum[£)[faraday2[[s, o]], coord[[s]]], {s, 1, 4}]+ Sum[aflBne[[i, i, j]]faraday2[|j, o]], 
{i, 1, 4}, {j, 1, 4}]+ Sum[affine[[o, fc, Z]]faxa<iay2[[A:, /]], {fc, 1, 4}, {Z, 1, 4}], {o, 1, 4}]] 



MatrixForm[CDfaraday2] 



/ \ 



-n[i,r-]+r(M(0'i)[t,r]+r(u(0'2)[t,r]-«(2,0)[(^r])) 



(A.12) 



This matches Equation (5.40). 



Now we must calculate (4V^ (FF^*^) + 7V^ (GF*^^)), which we shall call "righthand- 
side." To this end, we define FF^''' and GF*^''' as variables "Ffaraday2" and "Gmaxwell2" 
respectively, then we take their covariant derivatives "CDFfaraday2" and "CDGmaxwell2." 
Finally, we combine these so as to calculate "righthandside." 
We define FF^^ by entering: 

Ffaraday2:=Ffaraday2 = Fspecialcase faraday2 

Taking the covariant derivative V ^{FF^'^): 
CDFfaraday2:= CDFfaraday2 = FuUSimplifyf 

Table[Sum[£)[Ffaraday2[[s,o]],coord[[s]]], {s,l,4}]+ Sum[affine[[i,i,i]]Ffaraday2[[j,o]], 
{i, 1, 4}, {i, 1, 4}]+ Sum[affine[[o, fc, /]]Ffaraday2[[fc, /]], {k, 1, 4}, {/, 1, 4}], {o, 1, 4}]] 

Defining GF*^"'■. 

Gmaxwell2:=Gmaxwell2 = Gspecialcase maxwell2 

Taking the covariant derivative {GF*^'^): 
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CDGmaxweU2:= CDGmaxwell2 = FullSimpUfy[ 

Table[Sum[£)[Gmaxwell2[[s, o]], coord[[s]]], {s, 1, 4}]+ Sum[aflane[[i, i, j]]Gmaxwell2[[7, o]], 
{i, 1, 4}, {j, 1, 4}]+ Sum[affine[[o, fc, Z]]Gma^ell2[[fc, /]], {k, 1, 4}, {/, 1, 4}], {o, 1, 4}]] 

We compute j (4V/, {FFf"") + 7V^ (GF*'^'^)) by entering: 

righthandside:= righthandside = FuUSimplify [^CDFfaraday2 + ^CDGmaxwell2] 



We note that Equation (5.41) checks out since: 
righthandside [ [1]] 





and: 



righthandside [ [2] ] 



moreover U is given by: 
righthandside[[3]] * (45rV (4^2)) 

-6u[t,rf + 2ru[t,r]^ {-3u'^^''^\t,r] + 3ru'^^'^\t,r] - ru^^'^\t,r]) + 
r\[t,r] (6M(°'i)[t,r]2 - 2v^'^'\t,r]^ - 5t;(i'°)[t, r]^- 
2m(1'0) [t, r] [t, r] + Aru^^'^^ [t, r] 

3rn(i'°) [t, r]i;(i'i) [t, r]) + ri;(°'i) [t, r]^ (2u(°'2) [t, r] + Sw/^'O) [t, r] 

M(°'i)[t,r]2(6ra(0'2)[t,r]-2™(2'0)[t,r] 
M(0'i)[t,r] (2t;(°'i)[t,r]2-6M(i'°)[t,r 
8ra(^'°)[t,r]M(^'^)[t,r] +t;(^'°)[t,r] (5t;(i'°)[t, r] + 3rt;(^'i)[t, r] 

ri;(°'i)[t,r] (4t;(0'2)[t, r] - 7t;(2'0)[t, r]) 
r (-7t;(°'2) [t, r]u(^'°) [t, r]t;(i'°) [t, r] + 
M(°'2)[t,r] (-2M(i'°)[t,r]2 + 5t;(i'°)[t,r]2) + 
2 (3m(1'°) [t, r]2 + [t, r]^) m^^.o) ^] + 4^(1,0) ^]^{i,o) ^]^{2,o) ^ 

and V is given by: 



+ 
+ 
+ 
+ 
+ 

+ 
+ 
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righthandside[[4]] * (45r^/ (4q;^)) 



u[t, r]2 (-2t;(°'^) [t, r] + 2ri;(°'=^) [t, r] + SrT;^^'^) [t, r] 

14^^(1,0) [^^ ^] ^ ^^(0,1) [^^ ^] (4^(0,2) ^] _ 7^(2,0) ^] 

2M(°'i)[t,r] (2t;(°'i)[t,r]+ 2rt;(°'2)[t, r] + 5rt;(''°)[t, r]) 

M(°'i)[t,r] (3rt;(i'°)[t,r]u(i'i)[t,r]- 2M(i'°)[t,r] {2v^^''^^[t,r] + 7rv^'^^^[t,r] 

ri;(°'i)[t,r] (4M(°'^)[t, r] - 7M(^'°)[t, r]) 
t;(°'i)[t,r]2 (6rt;(°'')[t,r] -2rt;(2.o)[t,r] 
[t, r]2 (6t;(°'^) [t, r] + 2rt;(°'2) [t, r] + Srt;^^'^) [t, r] 
r (i;(°'2)[t,r] (5n(i'°)[t, r]^ - 2t;(i'°)[t, r]^ 

^(1,0) [^^ ^]^(l,0) ^] (_7^(0,2) [^^ ^] ^ 4^(2,0) ^] 



+ 

+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



2 (M(i'°)[t,r]2 + 3t;(i'°)[t,r]2) t;(2.o)[t, r]) (A.13) 



A. 5 The Maxwellian approximation 



Our next task is to check Equations (5.55) - (5.58). Introducing constants U and V (not to 
be confused with the functions U and V given above), we have that an eUiptically polarized 
ingoing cyhndrical wave is given by M[t,r] and f [t, r], where: 



u[t,r]:=^(BesseU[l,u}r]Cos[ujt\ - BesselY[l,a;r]Sin[a;t]) 



i>[t,r]:= ^(BesselJ[0,a;r]Cos[a;t]— BesselY[0,a;r]Sin[wt]) 

Note that the function u[t,r] now reverts back to denoting u (= Aq) again, instead of u. 
In order to calculate dr/dt for ingoing radial null geodesies, the equations of Section 
require us to calculate {{dtu — druY + r'^{dtv — drvY). Thereby we enter the lines: 



5.1.3 



FullSimpUfyf 

^ {{D[u%rlt] - D[u%rlr]f+ r\D[v%rlt] - D[v%Tlr]f)] 
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[/^((BesselJ[0,ra;] + BesselY[l, rcj])CosH + 
(BesselJ[l,rw] - BesselY[0, ra;])Sin[tu;])2 + 
V\{-BesseU[l,ru] + BesselY[0, ru;])CosH + 

(BesselJ[0, rco] + BesselY[l, ruj])Sm[tuj]f (A.14) 



This output (A.14) confirms Equation (5.55). 

For the outgoing radial null geodesies, we need ^ {{dtu + dru)'^ + r^{dtv + drv)'^). To 
this end, we enter: 

PiillSiinpUfy[ 

^ {iD[u[t,r],t] + D[u[t,r],r])H r^iD[v[t,r],t] + D[v[t,r],r]f)] 



f/^((-BesselJ[0, ru] + BesselY[l, ruj])Cos[tLu] + 
(BesselJ[l,rw] + BesselY[0, ra;])Sin[ta;])2 + 
V^{{BesseU[l, ru] + BesselY[0, ru;])Cos[tw] + 

(BesselJ[0, rcu] - BesselY[l, ruj])Sm[tu]f (A.15) 



This confirms Equation (5.56) 



Henceforth, we restrict to the case of circular polarization, whereby U = V = A = constant. 

U:=A 

V:=A 

We check Equations (5.57) and (5.58) by re-entering: 
FullSimpUfy[ 

^ {{DMt.rU - D[u[t,r],r]f+ r-'iDivit^rU - D[v[t,rlr])^)] 

A^i — + BesselJ[0, ruf + BesselJ[l, ruj?+ 

\ nrco 

BesselY[0,ra;]^ + BesselY[l, rw]^ ) (A.16) 



FullSimpUfyf 

^{{D[u[t,r],t]+D[u[t,r],r]f+ r\D[v[t,r],t] + D[v[t,r],r]y)] 
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A' 



+ BesselJ[0, rcof + BesselJ[l, rujY+ 



BesselY[0, ru]^ + BesselY[l, ru)\ 



(A.17) 



The outputs (A. 16) and (A.17) confirm Equations (5.57) and (5.58). 
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Appendix B 

Some additional remarks 



In Section 5.3.1 , we studied effective geometries corresponding to branch I of our exact 



solution B + 8a'^{E'^B — B^)/45 = k/r (5.71), witli certain restrictions on tlie constant E. 
The purpose of this appendix is to take a further look at the effective geometries of our exact 
solution. We do so only briefly. This appendix is intended to be read as a continuation of 



Section 5.3.1; we still assume that k > 0, and the (+) and (— ) polarization states have the 



same meaning as in Section |5.3.1 

Theorem 11. For E = 0, in the effective geometry of branch I corresponding to the (+) po- 
larization state, the outgoing radial null geodesies, issuing from any point where the effective 
geometry is defined, are never trapped. 

Proof. For the (+) polarization state, with E 



1 

A 



495 + 2Aa^B^ + 



0, Equation (5.74|) gives: 
135 - 136a^B^ 



224^2 



(B.l) 



We note that there is a particular radius vi such that for r > ri, we have B^ < 135/ [ISGa"^), 



and for < r < ri, we have 135/(136q;2 
have Bs = 15 /{8a^). 



<B' < Bi. Since E 



So for r > ri Equation (B.l) gives: 

1 

A r>ri 

Note that 1/A|,.>^^ > 0. 



45 
16a2 



-B' 



in the present case, we 



(B.2) 



Using Equation (5.80), we get that: 



dr 
di 



-EB 




r>r] 



r>ri 



(52 - E^) 



out, r>ri 



i?2 + l 



r>ri 



1-B^A 






r>ri 



(B.3) 
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For r > ri, the value of never exceeds 135/(136«^) (< B^), and we observe that this 
fact imphes dr/dt = ^jl — B'^A > 0. 



out, r>ri 



r>ri 



In the region where < r < ri, we have: 



1 

A 



45 5 
+ -5^ 



s<r<ri 28a^ 7 



So: 



di 



out, rs<r<ri 



1 



rs<r<ri 

For Ts < r < ri, the value of B^ never exceeds B^, and consequently dr/dt 
> 0. 



(B.4) 



(B.5) 



out, rs<r<ri 



1-B^A 



Ts<T<r\ 



□ 



Theorem 12. For E"^ < -^^^j in the effective geometry of branch I corresponding to the 
(— ) polarization state, either dr/dt out, r=rs is zero (when E > 0) or dr/dtm, r=rs is zero (if 
E < 0). If E = 0, then both dr/dt out, r=rs <^nd dr/dtin, r=rs o'^e zero. 



Proof. Since , we have that for the (— ) polarization state, at r = r^: 



1 

A 



Using (5.80) and (B.6), we get: 



and: 



dr 
di 



dr 
lit 



out, r = Ts 



= Bl 

B,{\E\~E) 
E^ + Bl ' 

B,{\E\+E) 



in, r = r„ 



E^ + B^ 



(B.6) 



(B.7) 



(B.8) 



If E > 0, then dr/dt 



out, r = Ts 



0, and dr/dt 



< 0. On the other hand, if E < 0, then 



m, r = r. 



dr/dt 
0. 



out, r = rs 



> 0, and dr/dt 



0. HE = 0, then dr/dt 



m, r = Te 



out, r = rs 



dr/dt 



□ 



Now let us say something about branch II. In this branch, B is an increasing function 
of r; it starts with the value of i?^ at r = r^, and increases towards the asymptotic value 
V^Bs at r = cxD. 

Theorem 13. For E ^ 0, in branch II, at infinity, the radial null geodesies propagate in 
only one direction; towards the axis if E > 0, away from the axis if E < 0. 
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Proof. In the limit r — )■ oo, we have B = \/^Bs, and Equation (5.74) gives: 



1 
A 



630 T 630 
224^2 

45 

or ^— ^ (depending on the polarization). 



(B.9) 



Equation (5.80) gives, for the radial null geodesies (both "ingoing" and "outgoing" 



dr 
dt 



_1_ 

A2 



45 
8a2 



E 



or 



45 
Aa2 



^2 + 



45 
8a2 



(depending on the polarization). (B.IO) 



Both polarization states travel inwards towards the axis if > 0; outwards to r = oo if 
E < 0. Note that this is consistent with Conjecture [2j We also note that, at infinity, 
the coordinate speeds of the two polarization states are reciprocal to one another and one 
polarization state propagates superluminally. 

□ 

In branch III, the absolute value \B\ is an decreasing function of r; at r = we have 
\B\ — )■ oo, and as r — )■ oo, we have \B\ — )■ y/SBg. At r = oo, the effective geometry 
corresponding to branch III is similar the effective geometry of branch II. 
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